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Abstract. Quantum lattice systems are rigorously studied at low temperatures. When the Hamiltonian 
of the system consists of a potential (diagonal) term and a — small — off-diagonal matrix containing 
typically quantum effects, such as a hopping matrix, we show that the latter creates an effective inter- 
action between the particles. 

In the case that the potential matrix has infinitely many degenerate ground states, some of them 
may be stabilized by the effective potential. The low temperature phase diagram is thus a small 
deformation of the zero temperature phase diagram of the diagonal potential and the effective potential. 
As illustrations we discuss the asymmetric Hubbard model and the hard-core Bose-Hubbard model. 
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1. Introduction 

Physics of a large number of quantum particles at equilibrium is very interesting and dif- 
ficult at the same time. Interesting, because it is treating such macroscopic phenomena as 
magnetization, crystallisation, superfluidity or superconductivity. And difficult, because their 
study has to combine Quantum Mechanics and Statistical Physics. 

A natural approach is to decrease difficulties arising from this combination by starting from 
only one aspect. Thus one can use only Quantum Mechanics and treat the particles first as 
independent, trying next to add small interactions. In the present paper we are concerned with 
the other approach. Namely, to start with a model treated by Classical Statistical Physics, 
adding next a small quantum perturbation. Another simplification is to consider lattice systems 
(going back to a physical justification for the modeling process, we can invoke applications to 
condensed matter physics). 

Quantum systems studied here have Hamiltonians consisting of two terms. The first term 
is a classical interaction between particles; formally, this operator is "function" of the position 
operators of the particles and it is diagonal with respect to the corresponding basis in occupation 
numbers. The second term is an off-diagonal operator that we suppose to be small with respect 
to the interaction. Typical example for this is a hopping matrix. 

The aim of the paper is to show that a new effective interaction appears that is due to 
the combination of the potential and the kinetic term. An explicit formula is computed, and 
sufficient conditions are given in order that the low temperature behaviour is controlled by 
the sum of the original diagonal interaction and the effective potential. To be more precise, 
it is rigorously shown that the phase diagram of the original quantum model is only a small 
perturbation of the phase diagram of a classical lattice model with the effective interaction. 

Thus, we will start by recalling some standard ideas of Classical Statistical Mechanics of 
lattice systems. The Peierls argument for proving the occurrence of a first order phase transition 



in the Ising model [Pei, Dob , Gri] marks the beginning of the perturbative studies of the low 
temperature regimes of classical statistical models. Partition functions and expectation values 
of observables may be expanded with respect to the excitations on top of the ground states, 
interpreting the excitations in geometric terms as contours. These ideas and methods are 



referred to as the Pirogov-Sinai theory; they were first introduced in [PS, 3in| and later further 



extended | Zah , BI, BS] 



The intuitive picture is that a low temperature phase is essentially a ground state config- 
uration with small excitations. A phase is stable whenever it is improbable to install a large 
domain with another phase inside. For such an insertion one has to pay on its boundary, it is 
excited (two phases are separated by excitations), but, on the other side, one may gain on its 
volume if its metastable free energy (its ground energy minus the contribution of small ther- 
mal fluctuations) is smaller than the one of the external phase. It is important to take into 
account the fluctuations since they can play a role in determining which phase is dominant. 
A standard example here is the Blume-Capel model with an external field slightly favouring 
the phase; at low temperatures, the "0" phase may be still selected because it has more 

low energy excitations (theory of such dominant states chosen by thermal fluctuations may be 



found in |BS]). 

The partition function of a quantum system Tr e~@ H may be expressed using Duhamel 
expansion (or Trotter formula), yielding a classical contour model in a space with one more 
(continuous) dimension. If the corresponding classical model (the diagonal part only) has 
stable low temperature phases, and if the off-diagonal terms of the Hamiltonian are small, the 
contours have low probability of occurrence and it is possible to extend the Peierls argument to 



quantum models | Gin | . More generally, one can formulate a "Quantum Pirogov-Sinai theory" 



[BKU, DFF1], in order to establish that (i) low temperature phases are very close to ground 



states of the diagonal interaction (more precisely: the density matrix A e @ H is close to the 
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projection operator \g)(g\, where \g) is the ground state of the diagonal interaction only) 
and (ii) low temperature phase diagrams are small deformations of zero temperature phase 
diagrams of the interactions. 

So far we have only discussed the case when the effect of the quantum perturbation is small, 
and the features of the phases are due to the classical interaction between the particles. It 
may happen, however, that the classical interaction alone is not sufficient to choose the low 
temperature behaviour. This is the case in the two models we introduce now and use later for 
illustration of our general approach. 

• The asymmetric Hubbard model. It describes hopping spin | particles on a lattice 
A C Z u . A basis of its Hilbert space is indexed by classical configurations n E {0, |, j, 2} A , 
and the Hamiltonian 

H = - ^cLcy^ + U^n^n^- fiJ2(n x ^ + n xl ) (1.1) 

||x-2/|| 2 =l er=T,j x x 

(the hopping parameter t a depends on the spin of the particle). In the atomic limit 
t+ = 1 1 = the ground states are all the configurations with exactly one particle at each 
site. The degeneracy equals 2' A ', which means that it has nonvanishing residual entropy 
at zero temperature. 

• The hard-core Bose-Hubbard model. We consider bosons moving on a lattice A C 
1? . They interact through an infinite on-site repulsive potential (hard-core), nearest 
neighbour and next nearest neighbour repulsive potentials. A basis of its Hilbert space is 
the set of all configurations n £ {0, 1} A , and its Hamiltonian: 

H = -t a x a y + Ui n x n y + U 2 n x n y - /i^n x . (1.2) 

||a:-J/|[a=l [|as-J/||2=l \\x-y\\ 2 =V2 x 

For U\ 2U2, and if < \i < 8U2, the ground states of the potential part are those 
generated by ( ) > i- e - an y configuration with each two lines empty, and the other anti- 
ferromagnetic, is a ground state (and similarly in the other direction). The degeneracy is 

of the order (if A is a square), there is no residual entropy. 

In these two situations, the smallest quantum fluctuations yield an effective interaction, and 
this interaction stabilizes phases displaying long-range order (there is neither superfluidity nor 
superconductivity) . 

In the case where classical and quantum particles are mixed in one model, like the Falicov- 



Kimball model, a method using Peierls argument was proposed by Kennedy and Lieb [KL|; it 



was extended in [LM] to situations that are not covered by the present paper, namely to cases 
of such mixed systems with continuous classical variables. 

Results very similar to ours have already been obtained by Datta, Fernandez, Frohlich 



and Rey-Bellet [DFFR]. Their approach is different, however. Starting from a Hamiltonian 
H(X) = if ( ) + XV, being a diagonal operator with infinitely many ground states, and V 
the quantum perturbation, the idea is to choose an antisymmetric matrix S = XS^ + X 2 S^ in 
such a way that the operator H^ 2 \X) = e s H(X) e~ s , expanded with the help of Lie-Schwinger 
series, turns out to be diagonal, up to terms of order A 3 or higher. If the diagonal part of 
if ( 2 ) has a finite number of ground states and the excitations cost strictly positive energy, it 
can be shown that the ground states are stable. It is possible to include higher orders in this 



perturbation scheme (see [ DFFR| ) 



In fact, our first intention was to study the stability of the results of [BS] with respect to a 
quantum perturbation, and we began the present study as a warm-up and the first simple step 
towards this goal. This simple step turned out however to be rather involved. Even though, 



at the end, the paper contains results similar to that of |DFFR], we think that the subject is 
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important enough to justify an alternative approach, and that there are some advantages in 
an explicit formula for the effective potential and sufficient conditions for it to control the low 
temperature behaviour that may be useful in explicit applications. 

The intuitive background of this paper owes much to the work of Bricmont and Slawny [BS] 



discussing the situation with infinite degeneracy of ground states, where only a finite number 
of ground states is dominating as a result of thermal fluctuations, and to the paper of Messager 
and Miracle-Sole ]MM| ] which was useful to understand the structure of the quantum fluctu- 



ations. Having expanded the partition function Tr e~@ H using Duhamel formula and having 
defined quantum contours as excitations with respect to a well chosen classical configuration, 
we identify the smallest quantum contours (that we call loops). Given a set of big quantum 
contours, we can replace the sum over sets of loops by an effective interaction acting on the 
quantum configurations without loops. This effective interaction is long-range, but decays ex- 
ponentially quickly with respect to the distance. This allows, for a class of models, to have an 
explicit control on the approximation given by effective interaction allowing to prove rigorous 
statements about the behaviour of original quantum model. 

An important model that does not fall into the class of models we can treat, is the (symmet- 
ric) Hubbard model. Take U = 1 and if = t \ = t in (|l.l|). Computing the effective potential 
stemming from one transition of a particle to a neighbouring site and back, we find an antifer- 
romagnetic interaction of strength t 2 . On the other hand, it is possible to make two transitions 
as a result of which the spins of nearest neighbours are interchanged, 

\n x ,n y ) = | |,|) = -c^Cy^c^ | t,|). 

It turns out that this brings the factor i 2 , which is of the same order as the strength of the 
effective interaction. In this case we cannot ensure the stability of the phases selected by the 
effective potential — we would need a stronger effective interaction. Otherwise the system 
jumps easily from a configuration with one particle per site to another such configuration, 
i.e. from a classical ground state to another classical ground state. We call quantum instability 
this property of the system. In the Hubbard model it is a manifestation of a continuous 
symmetry of the system, namely the rotation invariance. 

In Section ^ the ideas discussed above are introduced with the precise definitions. The 
effective potential is written down in Sections |2.2| (a general formula) an d|2.3| (a simpler formula 
in special cases). The results of the paper are summarized in Theorems |2.2| (a characterization 
of stable pure phases) and |2.3| (the structure of the phase diagram); experts will recognize 
standard formulations of Pirogov-Sinai theory. Taking into account that our aim is to describe 
in rigorous way the behaviour of a quantum system, some care must be given to the introduction 
of stable phases. We define them with the help of an external field perturbation of the state 
constructed with periodic boundary conditions. In Section |3| we apply the results to our 
two illustrative examples. The rest of the paper is devoted to the construction of a contour 
representation (Section ||), the proof of the exponential decay of the weights of the contours 
(Section ||), and, finally, the proofs of our claims with the help of contour expansions of the 
expectation values of local observables and the standard Pirogov-Sinai theory (Section ||). 

Let us end this introduction by noting that given a model which enters our setting, it is not a 
straightforward task to apply our theorems. One still has to separate the correct leading orders 
that determine the behaviour of effective interaction. This situation has the utmost advantage 
that it should bring much more pleasure to users, since the most interesting part of the job 
remains to be done — to get intuition and to understand how the system behaves. 

Acknowledgments. We are thankful to Christian Gruber for discussions. R.K. acknowl- 
edges the Institut de Physique Theorique at EPFL, and D.U. the Center for Theoretical Study 
at Charles University for hospitality. 
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2. Assumptions and statements 

2.1. Classical Hamiltonian with quantum perturbation. Let Z", v > 2, be the hyper- 
cubic lattice. We use \x — y\ '. — \\x — y||oo to denote the distance between two sites x, y G Z^. 
0, is the finite state space of the system at site x = 0, |0| = <S < oo. Our standard setting will 
be to consider the system on a finite torus A = (Z/LZ) 1 ' (i.e. a finite hypercube with periodic 
boundary conditions). With a slight abuse of notation we identify A with a subset of 7L V and 
always assume that it is sufficiently large (to surpass the range of considered finite range in- 
teractions). A classical configuration ra A (occasionally we suppress the index and denote it n) 
is an element of If A C A, the restriction of n A to A is also denoted by ua- 'Hk is the 
(finite-dimensional) Hilbert space spanned by the classical configurations, i.e. the set of vectors 

\v) = ^«n A l n A), a nA G C, 

with the scalar product 



"A 



( v \v') = ^2a* nA a' nA . 

n A 

Given two configurations ua G £l A and n' A , G £l A> , with A Pi A' = 0, it is convenient to define 
nAn' A , E Q AuA to be the configuration coinciding with ha on A and with n' A , on A' . 

The Hamiltonian is a sum of two terms, H F er = H ? per + V^ 01 . The former is the quantum 
equivalent of a classical interaction, the latter is the quantum perturbation. 

Assumption 1. Classical Hamiltonian. 

There exists a periodic interaction $ (i.e. a collection of functions <&a '■ Q A — > of finite 
range Ro (i.e. <&a = whenever diam^l > Rq) and period t§ such that 

< )per |n A ) = < )pcr (n A ) |n A ) = £ * A (nx) |n A >; 

AcA 

/or any foras AcZ" o/ side L i/iai is a multiple of £q and any n A 6 f2 A . 

When stressing the dependence on the interaction we will also use the notation F A 0)por (n A ) 
//'('""'("A). 

Let us suppose that a fixed collection of reference configurations G C is given[]. For any 
n G VtF" , ifZ", and any (finite range, but not necessarily translation invariant) interaction 
Vl/, we use e~(n) to denote the local contribution to the "energy" ^ of the configuration n at 
the site x, e*(n) = YIa^x ^ A \a\ A ^ ■ Notice that i/* per (n A ) = X^eA e z( n ) ^ or ever y n £ ^ A 
and A sufficiently large. Finally, let U{x) = {y G Z^; |y — x\ ^ -Ro}> ^4 = U X £aU(x) and 

We assume that the local energy gap of excitations is uniformly bounded from below, while 
the spread of local energies of reference states is not too big: 

Assumption 2. Energy gap for classical excitations. 

There exist strictly positive constants Aq and do such that: 

• For any x G Z u and any £ G(jr x ) , one has the lower bound 

e*(n) - maxe*(o) ^ A , (2.1) 



1 In some situations G is simply the set of all ground configurations of $. When discussing the full phase 
diagram, however, we will typically extend the interaction $ to a class of interactions by adding certain "external 
fields" . The set G then will actually play the role of ground states of the interaction with a particular values of 
external fields (the point of maximal coexistence of ground state phase diagram). 
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• and, 

m^\et(g)-et(g')\^8 . (2.2) 

Furthermore, we assume the following extension property on the set of reference states G: if, 
for a connected A C Z", a configuration n is such that njjf x \ G Gum for any x G A, then 
n A G G A . 

For later purpose, we note the following consequence of this assumption. 

Property. Let <E> satisfy Assumption R be such that R v ^ Aq/5q, and A <Z "L v with 
diamA ^ R. Then any pair of configurations g A G G A and n A £ G A , with n A \ A = g A \A' 
satisfies the lower bound 



£ 

A'cA 



[*A'(«aO " $A>(9A>)] > R~ v &a- (2.3) 



Proof. We choose g 1 G G such that g' A = g A . Then we have 

£ [t> A >(n A ,)-<f> A ,(g A ,)] =^(et(n A g' A c)-et(g')). (2.4) 

A'CA xCA 

Since n A £ G A , there exists at least one site x £ A such that njjr x -\ £ Gyi x \. From the assumption, this implies 
that 

^2 \$ A ,(n A >) -&A>(gA>j\ > A - 5 °- 

A'CA y£A,y^x 

Using \A\ ^ R v , we obtain the property. □ 

The quantum perturbation Vj^ eT is supposed to be a periodic quantum interaction. Namely, 
Vj^ er is a sum of local operators Va, V^ er = a Va, where Va has support supp A = A C A 
and A is, in general, a pair (A, a), where the index a specifies Va from a possible finite set 
of operators with the same support. We found it useful to label quantum interactions Va not 
only by the interaction domain A, but also, say, by quantum numbers of participating creation 
and annihilation operators. Thus, for example, the term A might, in the case of the Hubbard 
model, be a pair (<x, y >, f) corresponding to the operator Va = c x *Cy^. We refer to A as a 
quantum transition. 

Assumption 3. Quantum Perturbations. 

The collection of operators Va is supposed to be periodic^, with period Iq, with respect to the 
translations of supp A. The interactions Va are assumed to satisfy the following condition, 
for fermions or bosons, respectively: 

• ( Fermions ) Va is a finite sum of even monomials in creation and annihilation operators 
of fermionic particles at a given site, i.e. 



V A - £ <■> '< 



Xi,ai ■ ■ ■ ^x k ,a k ^yi,a 1 ■ ■ ■ ^ye,cr e 
(xi,ai),...,(xk,ak) 

with Xi,yi G A and o~i, a\ are the internal degrees of freedom, such as spins; k + 1 must 
be an even number. The creation and annihilation operators satisfy the anticommutation 
relations 



2 By taking the least common multiple, we can always suppose the same periodicity for <E> and V. Moreover, 
whenever a torus A is considered, we suppose that its side is a multiple of to. 
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• (Spins or bosons) The matrix element 

(nA\V A \n' A ) 

is zero whenever n^\A 7^ n ^\A an ^ otherwise it depends on ha and n' A only. 

In both cases V is supposed to have an exponential decay with respect to its support: defining 
the norm of a (quantum) interaction V b^ 

ll V H = F^\ v m , ax A \( n A\ v A\riA)\ , (2.5) 

we assume that ||V|| < 00. 

When stating our theorems, we shall actually suppose ||V|| to be sufficiently small. Notice 
also that we do not assume that V is of finite range, the exponential decay suffices. 

2.2. The effective potential. It is actually a cumbersome task to write down a compact for- 
mula for the effective potential in the general case. A lot of notation has to be introduced, and 
one pays for the generality by the fact that the resulting formulas look rather obscure; never- 
theless, the logic behind the following definitions and equations appears rather naturally along 
the steps in Section |||. In the next subsection we shall discuss a special case where the effective 
interaction is due to at most four transitions resulting in much simpler and straightforward 
formulas. We would like to stress that for typical concrete models this is entirely sufficient. 
The reader might thus skip the present subsection on the first reading and consider only the 
simplified situation of the next subsection. 

The real meaning of the next definitions (in particular, Q2.7| )) will appear more clearly only 
in Section but, in general case, we cannot leave it aside. First of all, we assume that a list S 
of sequences of quantum transitions A is given to represent the leading quantum fluctuations. 
The particular choice of S depends on properties of the considered model. Often the obvious 
choice like "any sequence of transitions not surpassing a given order" is sufficient. In general 
case, certain conditions (specified later in Assumption 5) involving S are to be met. For any 
gA £ Ga, the effective potential ^ is defined to equal 



— Ea E E 



n^l k 1 ,...,k n >2 {A \,...,Al i ,Al...,A% n )eS 

U iyj A*=A 

\\\ I { A i,---, A k 1 ',nAg A \A,---,nf 1 ~ 1 g A \ A )[]J(n^~ 1 \V A ^ 



fei-1 



cc<t,'<...<t;.<°° 



drj .. .dr|.[ Y[ e-^- T ^A^A^A>in£)-3> A ,{9 A ,)] j j 



I Tmin,- rj < and max,- ^ rj > Ol „ 

^ s : ^ VB ls ...B n ). (2.6) 

maxjj rj - mnijj rj 

To begin to decode this formula, notice first that the second sum is over all sequences 

(A\, . . . , Ah , A\ , . . . , A% n ) of transitions that are in the list S and are just covering the set 

A, UijA l j = A. The sum in the braces (for a given i = 1, . . . ,n) is taken over collections of 

configurations n A , . . . , n A l ~ l £ Ga with n A = n A l = gA, while the integral is taken over 
"times" attributed to transitions, with the energy term in the exponent taken over the set 
A i = U k j t 1 A),A = U xeA U(x). 



3 This is a norm, provided the multiplication of an interaction V by a scalar A is defined to be (XV) a = A' A ' Va • 
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Finally, there are some restrictions on the sums and integrals encoded in functions 
I [mini,j r)<0 and maxij rj>o] T . . 41 t,Jfei-l \ T , 

J„ r -- m ,,^ » V ' ' ' B ")> and ■ • • ' A fe 4 ; n i 5A\A, • • • , nX 9a\a)- The 

easiest is the first one. One just assumes that the interval between the first and the last of 
concerned "times" contains the origin and the integrand is divided by the length of this interval. 
The function ^(Bi, . . . B n ) in terms of the sets Bi = A 1 x [r{, r^J cZ"x [—00, 00], i = 1, . . . ,n, 
is the standard factor from the theory of cluster expansions defined as 



<p T (Bi,...,B n ) 



1 if n = 1 

Eg Ue(ij)eg (~ 1 [ R i u B J is connected ] ) if n ^ 2 



with the sum over all connected graphs Q of n vertices. Connectedness of a set B d 7L V y. 
[—00,00] is defined by combining connection in continuous direction with connection in slices 
{x\(x,t) G B} C 7L v through pairs of sites of distance one. The most difficult to define is the 
restriction given by the function X that characterizes whether the collection of transitions is 
connected, in some generalized sense, through the intertwining configurations. A consolation 
might be that in lowest orders it is always true. Namely, I(A\, . . . , A).; n\, . . . , n 1 ^ 1 ) = 1 
whenever k ^ 5. To define it in a general case, consider A\, . . . y A^ C 7H 1 and n 1 , . . . , n fc_1 G 
VtF . Taking A = U x& aU(x) and E{n) = {x G A : njjr x \ ^ gjj{x) f° r an Y 9 € G}, we consider 
the set C Z" +1 , 



fe r 



#°> = U 
3=1 



fc-1 



Aj x {2j - 2} U LJ i £(n J ) x {2j - 1} 



Think of layers, one on top of another — configurations on odd levels interspersed with transi- 
tions on even levels. The set B^ decomposes into connected components, B^ = Lie^i B^°\ 
To any B^\ define the box B^ C W +l as the smallest rectangle containing B^\ Then 
let B^ = Lie ^ 1 Bjp ; decompose into connected components B^ = Lie 1 and repeat 

the procedure until no change occurs any more, i.e. until £?( m ) = Lie > 1 B^ m \ The function I 
characterizes whether this final set, the result of the above construction, is connected or not, 

~, . , 1 fr_is I 1 if i?( m ) is connected , 

l(A 1 ,...,A k ;n\,...,n k A 1 ) = { . (2.7) 

otherwise. 



2.3. Quantum fluctuations with less than four transitions. The equation ( [2.6D for the 
effective potential is hard to handle in general case. However, in many situations it is enough 
to consider only small sequences of less than four quantum transitions to define it. We rewrite 
in this section the explicit formula? for the effective potential in such a case. 

We assume thus that a list S of sequences of quantum transitions A, containing at most 4 
transitions, is given to represent the most important quantum fluctuations. Let us decompose 
S = u U S( 4 \ with denoting the list of sequences with exactly k transitions, and 
write 

_ i|,(2) _|_ ^(3) _|_ ^(4)_ ^.8) 

Here is the contribution to the effective potential due to the fluctuations from S^ k \ 
Let 

4>A{n A ]g A )= ^ {$A>(n A ,) - $a'(9a>) ■ 
A'cA 



s 
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Then, for any connected A C 1> u and qa £ Ga, we define 

E E ^S^ . <™> 

AiUA 2 =A 

\_ (gAj \ n A)(nA\ Va 2 \n' A )(n' A \ V Az \9A) , 0ln x 

(Ai,A2,A3_)eS (3) n A ,n^G A 
AiUA 2 UA 3 =A 



'M^a^a^a^^S'a) 



The expression for ^( 4 ) becomes more complicated (we shall see in Section || that clusters of 
excitations are actually occurring here), 



(4i,A 2 ,A3_,A4]e5( 4 ) 
AiUA 2 UA 3 uA 4 =A 



E 



(gAj \&! \n A )(n A \V A2 \n' A ){n' A \ V A ., \n' A )(n' A \ V Aa \g A ) 



<t>A ( n A ;9A)<t>A (n' A ;sa)0a (n A ;g A ) 
n A ,n' A ,n' A ^G A 



\ (3.41^! |n A )(n A |\^ 2 |g A )( 3A |V^ ? \n' A )(n' A \V A4 \g A ) f i i ^ 



2 

n A ,n^G A 



. (2.11) 



Above we denoted A 1 = A\ U A 2 and yl 2 = I3 U ^4. Property ( |2.3| ) implies that all the 
denominators are strictly positive. 

These equations simplify further if Va is a monomial in creation and annihilation operators; 
indeed in the sums over intermediate configurations only one element has to be taken into 
account. 

Notice, finally, that the diagonal terms in V are not playing any role in the previous defini- 
tions; we consider that they are small, since otherwise we would have included them into the 
diagonal potential. 

2.4. Stability of the dominant states. The aim of rewriting a class of quantum transitions 
in terms of the effective potential was to get a control over stable low temperature phases. To 
this end, the three conditions, expressed first only vaguely and then in precise terms in the 
following Assumptions 4, 5, and 6, must be met. Namely, we suppose that 

• the Hamiltonian corresponding to the sum <3? + of the classical (diagonal) and effective 
interactions has a finite number of ground configurations, and its excitations have strictly 
positive energy^]; 

• the list S contains all the lowest quantum fluctuations; 

• there is no "quantum instability"; the transition probability from a "ground state" g to 
another "ground state" g' is small compared to the energy cost of the excitations. 

Each component of the effective interaction *&a is a mapping Ga —> R; let us first extend 
it to £l A — > R by putting *a(^a) = if n A ^ Ga- To give a precise meaning to the first 
condition, we suppose that a finite number of periodic reference configurations D C G is given 
such that the interaction + ^ satisfies the Peierls condition with respect to D. We choose a 
formulation in which it is very easy to verify the condition and, in addition, it takes into account 
the fact that the configurations from D are not necessarily translation invariant. Namely, we 
will formulate the condition in terms of a potential T that is equivalent to $ + and is chosen 
in a suitable way. Of course, in many particular cases this is not necessary and the condition as 
stated below is valid directly for $ + ^. However, in several important cases treated in Section 
3, the interaction $ + turns out not to be so called m-potential and the use of the equivalent 



4 Again, when exploring a region of phase diagram at once, we have a fixed finite set of reference configurations 
that, strictly speaking, turn out to be ground configurations of the corresponding Hamiltonian for a particular 
value of "external fields". See below for a more detailed formulation. 
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m-potential T not only simplifies the formulation of the Peierls condition, but also makes the 
task of its verification much easier. 

We will consider the interactions tp and cj) to be equivalent if, for any finite torus A and any 
configuration n G A , one has 

Hl pc \n) = Ht per (n). 
Notice that the above amounts also to the equality 



y^enn) - > ( ,.(//) 

xeA ireA 



Assumption 4. Peierls condition. 

There exist a finite set of periodic configurations D C G with the smallest common period Lq, 
a constant A such that A > || V\\ k for some finite constant k, and a periodic interaction 
T = {T^} (with period £q) that is equivalent to <3> + ^ such that the following conditions are 
satisfied. The interaction T is of a finite rangef\ R G N such that 

R v ^ AqA), (2.12) 

with the constants 5o and Aq determined by the interaction <3? in Assumption |^. The value 
ej(c?) is supposed to be translation invariant with respect to x for any d G D, and the 
interaction T satisfies the following conditions: 

• For any x G A and any n with Hum ^ Gjjt x \, one has 

ej(n) - maxej(s) ^ ±A . 

• For any x G A and any n with ny^ ^ Dy^, V(x) = {y G A; \y — x\ < R}, one has 

e I( n ) -minej(d) ^ A. 

The following assumption is a condition demanding that the list S should contain all tran- 
sitions that are relevant for the effective potential. We define 

m(V Al ,...,V Ak ) = max max \(g\ V Al \n l )(n l \ V M |n 2 ) . . . (n k - l \ V Ak \g)\. (2.13) 

g&G n 1 ,...,n k - x $G 

Assumption 5. Completeness of the set of quantum transitions. 

There exists a function bi(-) with lim^o bi(\) = such that for any sequence 
(Ai, . . . , A m ) $l S with connected U™ ^Ai one has 

m(V Al , . . . , V Ak JmiV^ , . . . , V Ak J . . . m(V Akn _ i+1 , . . . , V Am ) <h(\\V\\)A. 

Finally, we have a condition assuring that there is no quantum instability. 

Assumption 6. Absence of quantum instability. 

There exists a function 62(0 with lim^o ^2 (A) = such that for any sequence (Ai, . . . , A m ), 
and any g,g' G G, g 7^ g' , one has 

(g\V Al ...V Am \g') 2 (||V1)A. 



5 Usual notion of (physically) equivalent interactions (see [Geo , |EFS | ) is slightly weaker, but we will not 
need it here. 

6 We will suppose, taking larger R if necessary, that it is larger or equal to the range Ro of $, as well as to 



the range of the effective interaction ^ and to Lq. 
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2.5. Characterization of stable phases. Notice first that the specific energy per lattice site 
of the configuration d € D, defined by 

e(d) = hm JL l$A(d A ) + *A(d A )}, (2.U) 
' ' AcA 

is equal, according to Assumption 4, to ej(d) (whose value does not depend on x). 

Our first result concerns the existence of the thermodynamic limit for the state under periodic 
boundary conditions. Taking Lq to be the smallest common period of periodic configurations 
from D, we always consider in the following the limit over tori A /* Z u whose sides are multiples 
of Lq and £q. 

Theorem 2.1. Thermodynamic limit. 

Suppose that the Assumptions^-^ are satisfied. There exist constants Eo > and (3q = (3q(A) 
( depending on u, S, R, £q ) such that the limit 



TrTe-P H 



A 



<T) pcr = hm —ftr- (2.15) 

'P A/Z" Tr e _/3ff A 

exists whenever \\V\\ ^ £o; P ^ Po> an d T is a local observable. 

Notice the logic of constants in the theorem above (as well as in the remaining two theorems 
stated below). We first choose eo- Then, for any ||V|| ^ eo one can choose (5q (depending 
on A that is determined in terms of V through the effective potential such that the claim 
is valid for the given V and any (3 ^ /3o(A). With ||V|| — > we may have to go to lower 
temperatures (higher 0) to keep the control. Of course, if A does not vanish with vanishing 



||V|| (i.e. Assumption 4 is valid for $ alone) as was the case in |BKU], one can choose the 
constant (3q uniformly in ||V||. 

If there are coexisting phases for a given temperature and Hamiltonian, the state (•) ( g Cr will 
actually turn out to be a linear combination of several pure states. A standard way how to select 
such a pure state is to consider a thermodynamic limit with a suitably chosen fixed boundary 
condition. In many situations to which the present theory should apply, this approach is not 
easy to implement. The classical part of the Hamiltonian might actually consist only of on-site 
terms and to make the system "feel" the boundary, the truly quantum terms must be used. One 
possibility is, of course, to couple the system with the boundary with the help of the effective 
potential. The problem here is, however, that since we are interested in a genuine quantum 
model, we would have to introduce the effective potential directly in the finite volume quantum 
state. Expanding this state, in a similar manner as it will be done in the next section, we would 
actually obtain a new, boundary dependent effective potential. One can imagine that it would 
be possible to cancel the respective terms by assuming that the boundary potential satisfies 
certain "renormalizing self-consistency conditions". However, the details of such an approach 
remain to be clarified. 

Here we have chosen another approach. Namely, we construct the pure states by limits of 
states (•)* per , defined by ( 2.15| ) with H^ er = per + V per , where <3? a is a perturbation of 



the interaction $ suitably chosen in such a way that one approaches the coexistence point from 
the one phase region. Consider thus Trq, the space of all periodic interactions of range i?o- 
We say that a state (•)« per , <f> S ^r , is unperturbable if it is insensitive to small perturbations: 



(T)^' per = lim(T}Jf +OT/,)per (2.16) 



for every ip £ and every local observable T. We define now a state (-)a to be a pure 
stable state (with classical potential $ and quantum interaction V) if there exists a function 
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(0, q?o) 3 oc — > 6 J r R() so that lim a _o+ 
for every local observable T. 



the states (•}? per are unperturbable, and 



a->0+ 'P 



(2.17) 



Theorem 2.2. Pure low temperature phases. Under the Assumptions [J^| and for any 

r) > 0, there exist eq > and /?o = Po(A) (depending on u, S, R, Iq) such that if \\V\\ ^ So and 
(3 ^ f3o, there exists for every d € D a function f^{d) such that the set Q = {d € D; Re /^(d) = 
mm d'e-D R e f^(d')} characterizes the set of pure stable phases. Namely, for any d € Q: 

a) The function f^(d) is equal to the free energy of the system, i.e. 

f(d) = -1 hm -j- log Tr e-^" . 
' w /3 A/Z" |A| 

b) There exists a pure stable state (•}^. Moreover, it is close to the state \d\) in the sense 
that for any bounded local observable T and any sufficiently large A, one has 



(T) d -(d A \T\d A ) 



^ 7/| suppT|||T| 



There is an exponential decay of correlations in the state {•)% 



where suppT is the support of the operator T. 
There is an expo 
£ d > such that 

(TT% - {T)*(T% 



i.e. there exists a constant 



s$ | suppT|| suppT'|||r||||T'|| e - dist ( su PP T - su PP T ')/« d 



for any bounded local observables T and T' . 



d) The state (•) ( 3 >cr is a linear combination of the states (•}% 



d E Q, with equal weights, 



for each local observable T. 



2.6. Phase diagram. We now turn to the phase diagram at low temperatures. Let r be the 
number of dominant states, i.e. r = \D\. To be able to investigate the phase diagram, we 
suppose that r — 1 suitable "external fields" are added to the Hamiltonian H per . Or, in other 
words, we suppose that classical potential and quantum interaction V depend on a vector 
parameter fi = (/ii, . . . ,fJ, r -i) £ U, where U is an open set of R r_1 . The dependence should 
be such that the parameters fx remove the degeneracy on the set D of dominant states. One 
way how to formulate this condition is to assume a nonsingularity of the matrix of derivatives 



dm 



Assumption 7. The potentials <3? and V are differentiable with respect to fi and there exist a 
constant M < oo such that 

d 

—®A{n A ) <M 
dm 



sup 



and 



< M 



for all fj, €LlA. 

Further, there exists a point /jq £W such that 



(d) = e W) (d') for all d,d' £ D, 
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and the inverse of the matrix of derivatives 




has a uniform bound for all \i $LU. 

Notice that if for some d S D one has e^(d) = e^ := min^/go e^(d'), then, according to 
Peierls condition (Assumption 4), the configuration d is actually a ground state of T. Thus, 
the assumption above implies that the zero temperature phase diagram has a regular structure: 
there exists a point jiQ € IA where all energies e fl0 (d) are equal, e fl0 (d) = e^°, r lines ending 
in hq with r — 1 ground states, ^r(r — 1) two-dimensional surfaces whose boundaries are the 
lines above with r — 2 ground states, . . . , r open (r — 1) -dimensional domains with only one 
ground state. Denoting the (r — \Q\) -dimensional manifolds corresponding to the coexistence 
of a given set Q C D of ground states by 

Wl*(Q) = \a G W;Ree^(d) = minRee'V) for all d € q), (2.18) 
I d'eD ) 

we can summarize the above structure by saying that the collection V* = {97T (Q)}qcd deter- 
mines a regular phase diagram. Notice, in particular, that Uq c d97T (Q) = U, TT (Q)nTT (Q') = 
whenever Q ^ Q' , while for the closures, Wf(Q) nM*(Q') = 971* (Q U Q'). Here we set 
971(0) = 0. 

The statement of the following theorem is that the similar collection V = {^01(Q)}qcd of 
manifolds corresponding to existence of corresponding stable pure phases for the full model is 
also a regular phase diagram and differs only slightly from V* . To measure the distance of two 
manifolds 971 and 971', we introduce the Hausdorff distance 

dist H (971,971') = max(sup dist (fi, 971'), sup dist (//, 971)). 



Theorem 2.3. Low temperature phase diagram 

Under the Assumptions [J-^ there exist £q > and Po = (3o(A) such that if 
ll^ll+Si=i llflr^ll ^ £ o and f3 ^ (3q, there exists a collection of manifolds V 13 = {97t^ (Q)}qcd 
such that 

(a) The collection determines a regular phase diagram; 

(b) If \x G 97l^(<5), the corresponding stable pure state (-)^ exists for every d € Q and satisfies 
the properties b), c), and d), from Theorem \2.Q ; 

(c) The Hausdorff distance dist h between the manifolds of and their correspondent in V* 
is bounded, 

r - 1 || ay 

dist H (97t$(Q),97^(Q)) <O(e-0 +\\V\\+J2 k~ )' 

for all Q C D. 

The proofs of these theorems are given in the rest of the paper. Expansions of the partition 
function and of expectation values of local observables are constructed, and interpreted as 
contours of a classical model in one additional dimension. Then we show that the assumptions 
for using the standard Pirogov-Sinai theory are fulfilled, and, with some special care to be 
taken due to our definition of stability, the validity of the three theorems follows. 



3. Examples 
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3.1. The asymmetric Hubbard model. The state space is = {0, f, |, 2} and the Hamil- 
tonian is written in ( |1 . 1|) . Hence the classical interaction is 

'0 ifn x = 

${x}(nx) = < -H if n x =| or n x =| (3.1) 
U-2fi iin x = 2 

(Rq = 0). We choose the chemical potential such that < fj, < U. The set G is here the set of 
ground states of <£, i.e. 

G = {ne VF : n x =] or n x =j for any x G Z v }. 

Assumption [2] holds with Ao = min(^, U — fi) and 5o = 0. 
The quantum perturbation is defined to be 

y A= Ul c U c v1 if A = (<x,y>,T) ^ 
and we always have A = {x,y} for a pair of nearest neighbours x,y G Z" . ||V|| = |ft|5 (if 

I*tI>I*iI)- 

The sequence S of transitions that we consider is 

S = {(A, A') : A = (< x, y >, f ) and A' = (<y, x>, f) for some x, y G Z 17 , ||x — 2/ 1 1 2 = 1}- 

The effective potential is given by Equation (|2.9|). For any x,y G Z", nearest neighbours, 
any configuration n such that \n) = c^cyp \g), g G G, has an increase of energy of 

</ , {x,j/}( ra { : r,y};9{a:,s/}) = ^- 

Furthermore we have 

•t.^, ,u/„, ,,,t, J, ^-J 1 if »K« } e{(T, l),(U)} 



(flfoy}! C a;t C S/T C j/t Cx T Ifffol/}) + (d{x,y}\ C y\ c x\C- x ^y\ \9{x,y}) 



otherwise. 



(3.3) 



Therefore 



Let V(x) = {y G Z^ : \y - x\ < 1}. Taking T = $ + we have 

||a^ — 2/|| 2=1 

(ej(ny( 2 .)) was initially a function — > R, but is actually depends on ny( x ) only). 

The set D has two elements, namely the two chessboard configurations S 1 ^ and d^; if 

(-ir := nr=i(-irs 

1 j, if (—1)* = — 1 x \i if c-i) x =i. 

The last inequality of Assumption [| holds with A = |t?/f7. 

The maximum of the expression in Assumption || is equal to 2U max(tf/tp tT). If there exists 
e > such that \t± | ^ |tf the bound of Assumption ||| can be chosen to be b\ = 2U\t^ \ 2e . For 
Assumption |6| the expression has maximum equals to 2Z7|ij|/|i|| and we can take 62 = 2C/|i-|-| £ 
(this Assumption is not true in the symmetric Hubbard model; the effective potential is not 
strong enough in order to forbid the model to jump from one g to another g'). 
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Our results for the asymmetric Hubbard model can be stated in the following theorem (first 
obtained by |pFF2fl ). 

Theorem 3.1 (Chessboard phases in asymmetric Hubbard model). Consider the lattice Z u , 
v ^ 2, and suppose < \i < U and \t\ \ ^ %| 1+e with e > 0. Then for any 6 > 0, there 
exist to > and /?o(%) < o° (hn^T^o A)(*f) = oo) such that if \t^\ ^ to and (3 *2 flo, 

• The free energy exists in the thermodynamic limit with periodic boundary conditions, as 
well as expectation values of observables. 

• There are two pure periodic phases, and (-}^\ with exponential decay of correlations. 

• One of these pure phases, is a small deformation of the chessboard state \d^}: 

w {>i-5 if(-ir = i s d)Us if(-ir = i 
{ xV p\^5 if (-ir = -i {xl) e \>i-s if(-ir = -i. 

The other pure phase, , is a small deformation of \d^}. 
To construct the two pure phases, one way is to consider the Hamiltonian 

Hr r (h) = Hr - h^(-lT(n x1 - n xi ). 

Then 
and 



where (•)l' er (/i) is defined by ( 2.15 ) with Hamiltonian iS"^ er (/i). 

3.2. The hard-core Bose-Hubbard model. The state space is f2 = {0, 1} and the Hamil- 
tonian is written in (|1.2|). Let P a plaquette of four sites; the classical interaction is 



®p(np) = \U\ n x n y + U 2 n x n y - i/x^n x , (3.6) 

\\x-y\\ 2 =l \\x-y\\ 2 =V2 

and = if A is not a plaquette. Here Ro = 1. Remark that we have 

$ P (np) = {\U X - \U 2 ) ^ + n v-\f + u ^{Y. n -- l 2-wS + C (3 ' 7) 

||a:-l/||a=l 

with C = —\U\ + ^U 2 — \n — /U 2 . $p(np) is minimum if np = ( J jj) , or an y configuration 
obtained from ( J {] ) by rotation. Hence we define 

G={nG{0,lf:upe{(ig),(8J),(8S),(»0)} for any plaquette p] 

(G is here the set of ground states of the interaction $). Since &p(np)—§p(gp) ^ \ min(//, &U 2 — 
fj,), for any np Gp, gp G Gp, Assumption [2] holds with Ao = jq min(//, 8U 2 — n), provided 
< 11 < 8U 2 . 8 = 0. 

We take as sequence of transitions for the smallest quantum fluctuations 

S = {(A, A') : A =<x, y > and A' =<y,x> for some x, y G Z 2 , \\x — y\\ 2 = 1}. 
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The effective potential follows from Let V xy = \Jp n i X y\^$P and more generally we 

denote by V any 3 x 4 or 4 x 3 rectangle. Up to rotations, we have to take into account five 
configurations, namely 

010 010 101 101 100 

ooo ooo ooo ooo 001 

010 101 010 101 100 

000 000 000 000 001 

(A) (B) (C) (-D) (E) 

9-p 9j> 9j> 9-p 9-p 

We find ^v{g { v ] ) = -t 2 /2U u ^v(gP) = ~t 2 /AU 2 , and Ms?) = = *^C#°) = 

0. 

The equivalent potential T should have a range bigger or equal to 3; in this case, the 
corresponding energies ej(n) would depend on the configuration on the square 7x7 centered 
at x. We could proceed in this way, but actually it simplifies a lot to do the following. Looking 
in the derivation of the contour model (Section 4.1), we see that beside of Assumption the 
only property that ej has to satisfy is that 



el(n) 



see equation ( 4.29 )]. Therefore we do not define T, but we do define 



el{n v{x) ) = \Y J $p{n P ) + \ £ W^) + t E (3-8) 

PBx y,\\y-x\\ 2 = l VCV(x) 

where V{x) is the square 5x5 centerd at x. 

The configurations gy(x) £ Gv(x) are ( U P to rotations and reflections) 

10101 10101 10101 10101 oioio ooooo ooooo ooooo 

00000 00000 00000 00000 00000 10101 10101 01010 

10101 10101 01010 01010 10101 ooooo ooooo ooooo 

00000 00000 00000 00000 00000 10101 01010 01010 

10101 01010 01010 10101 01010 ooooo ooooo ooooo 

W 00 (c) (d) (e) (/) (g) (h) 

9v(x) 9v(x) y V(x) 9y( x ) 9v(x) y V(x) »V(i) 

We find that ej(^ } ) = ej(^ x) ) = ej^) = -t 2 ^!"* 2 /^, otherwise ej( 3y(;c) ) ^ - 
3t 2 /4f/i -t 2 /8U 2 (if C7i ^ 2C/ 2 ). This allows to take A = t 2 (^ - ^-) in Assumption |. The set 
of dominant states D consists in all the configurations generated by gyl x \ and gy) x y \D\ =8. 

The maximum of the expression in Assumption || is b\ = i 2 (g^ ~~ i^")" 1 - m Assumption |6| 
b 2 = 0, because g ^ g' means that g and g' must differ on a whole row, and the matrix element 
is zero for any finite m. 

These eight dominant states bring eight pure periodic phases, , . . . , ; each one 
can be constructed by adding a suitable field in the Hamiltonian (e.g. the projector onto the 
dominant state). 

Theorem 3.2 (Hard-core Bose-Hubbard model). Consider the hard-core Bose-Hubbard model 
on the lattice 7?, and suppose U\ > 2U 2 and < /j, < 8U 2 . There exist to > and [3o(t) < oo 
dim^o Po(t) = oo) such that ift^to and (3 ^ (3q, 

• the free energy exists in the thermodynamic limit with periodic boundary conditions, as 
well as expectation values of observables, 

• there are 8 pure periodic phases with exponential decay of correlations. 



Each of these eight phases is a perturbation of a dominant state d, and the expectation value 
of any operator is close to its value in the state d, see Theorem |2.2j for more precise statement. 
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4. Contour representation of a quantum model 

Our Hamiltonian has periodicity Iq < oo. Without loss of generality, however, one can 
consider only translation invariant Hamiltonians, applying the standard trick. Namely, if Q is 
the single site phase space, we let O' = Qi^'—M". g' = |jy| = 5^0. Then we consider the torus 
A' C Z u , £q\A'\ = |A|, each point of which is representing a block of sites in A of size £q, and 
identify 

n ,A ' ~ n A . 

Constructing TL' as the Hilbert space spanned by the elements of /A , it is clear that TL' 
is isomorphic to TL. The new translation invariant interactions <£' and V are defined by 
resumming, for each A C A', the corresponding contributions with supports in the union of 
corresponding blocks. Notice the change in range of interactions. Namely, it decreased to 
\R/£o] (the lowest integer bigger or equal to R/£q). 

From now on, keeping the original notation TL, S, . . . , we suppose that the Hamiltonian is 
translation invariant. 

The partition function of a quantum model is a trace over a Hilbert space. But expanding 
e~@ H with the help of Duhamel formula we can reformulate it in terms of the partition function 
of a classical model in a space with one additional dimension (the extra dimension being 
continuous). In this section we present such an expansion, leading to a contour representation, 
of the partition function Z? ei := Tr e~P H A in a finite torus A. 

Expansion with the help of Duhamel formula yields 

e -W=£ £ f dn...dr m 
m > A u ...,A m J 0<T 1 <...<r m </3 
AiCA 

e- rifJ A V Al <r {T2 ~ Tl)H * V A2 ...V Am e~^- T ^ H ^ . (4.1) 
Inserting the expansion of unity 1% A = ]Cn A I^a) {^a I to the right of operators YAj-, we obtain 

m > n l ,... n m A lr ..,A m ^<ri<-<r m <P 
AiCA 

e -nHf-K) (n i , jn 2 } e -(^-n)flf . . . (< | ^ |4 } e -(/^ m )<>-«) . (4 . 2) 

This expansion can be interpreted as a classical partition function on the (y + l)-dimensional 
space A x [0, /?]. Namely, calling the additional dimension "time direction", the partition 
function Z^ CT is a (continuous) sum over all space-time configurations = ti\{t), t £ [0, /3], 
and all possible transitions at times corresponding to discontinuities of n\{r). Notice that 
72a (t) is periodic in the time direction. Thus, actually, we obtain a classical partition function 
on the d + 1-dimensional torus Ta = A x [0, f3] per with a circle [0, /3] per in time direction 
(for simplicity we omit in Ta a reference to (3). Introducing the quantum configuration wt a 
consisting of the space-time configuration ti\(t) and the transitions (Aj,Tj) at corresponding 
times, we can rewrite ( |4.2j ) in a compact form 

Zf" = J du TA p^(u TA ) (4.3) 



with p peT (uj A ) standing for the second line of (4.2). 



EFFECTIVE INTERACTIONS DUE TO QUANTUM FLUCTUATIONS 



17 



Now, we are going to specify excitations within a spacetime configuration n and identify 
classes of small excitations — the loop^\ — and large ones — the quantum contours. 

A configuration n G is said to be in the state g G G at site x whenever njjM = gjjix) 
(notice that, in general, g is not unique). If there is no such g £ G, the configuration n is said 
to be classically excited at x. We use E(n) to denote the set of all classically excited sites of 
n G For any AcF, let us consider the set of quantum configurations on the torus 
Ta- Whenever u G Q A , its boundary B^'{u) C Ta is defined as the union 

JJ<°>(«) = (U TeM (E(n(r)) x r)) U (U™^ x tj)). (4.4) 

The sets A\ x C Ta represent the effect of the operator V and for this reason are called 
quantum transitions. It is worth to notice that the set B^°\u) is closed. 

Next step is to identify the smallest quantum excitations — those consisting of a sequence of 
transitions from the list S. First, let us use B^°\uj) to denote the set of connected components 
of B (0) (cj) (so that B (0) (w) = U Be0 (o) (a ,)-B). To any B G £ (0) (cj) that is not wrapped around 
the cylinder (i.e., for which there exists time tb G [0, (3} per with B n (Z u x tb) = 0) we assign 
its sequence of transitions, S(B,u:), ordered according to their times (starting from tb to [3 
and proceeding from to tb) as well as the smallest box B containing B. Here, a box is any 
subset of Yjv of the form A x [ti, T2] with connected AdTf and [n, T2] C [0, /3] per (if t\ > T2, 
we interpret the segment [ti,T2] as that interval in [0, /3] per (with endpoints t\ and t<i) that 
contains the point = /3). 

We would like to declare the excitations with S(B, u) G S to be small. However, we need to 
be sure that there are no other excitations in their close neighbourhood. If this were the case, 
we would "glue" the neighbouring excitations together. This motivates the following iterative 
procedure. 

Given u, let us first consider the set B^\u) of those components B G B^(u) that are not 
wrapped around the cylinder and for which S(B,uj) G S, where S is the set of all subsequences 
of sequences from S. Next, we define the first extension of the boundary, 

= ( U B £B (0) M <) M 5 ) U ( U B £ <>(^)- 

Using £>W(cj) to denote the set of connected components of B^\u) and B^\u) C B^^u) 
the set of those components B in B^(u)) that are not wrapped around the cylinder and for 
which| S(B,u) G S, we define 

Iterating this procedure, it is clear that after finite number of steps we obtain the final 
extension of the boundary, 

B(u) = (U BgBW(wAB w (w) B) U (U BeB ( k)(w) S). 

Here, every B G B^ (cj) is actually a box of the form Ax[n, tq\ (that is not wrapped around the 

cylinder) and S(B,u) G S. Let us denote B(u) = B^\u) and consider the set i3o(w) C B(u>) 

of all those sets B G B^\u) for which actually S(B,u) G S and, moreover, ua{t\ — 0) = 
^a(t"2 + 0). Finally, let Bi(u) = B(lj) \ Bo(u) — it represents the set of all excitations of 
uj that are not loops. Taking, for any closed B C Ta, the restriction n# of a space-time 



7 Even though the present framework is more general, the name comes from thinking about simplest excitations 
in Hubbard type models. Namely, a jump of an electron to a neighbouring site and returning afterwards to its 
original position. 

8 A set B 6 B^\uj) may actually contain several original components from £>g° (oi). We take for S(B,lj) the 
sequence of all transitions in all those components. 



18 



EFFECTIVE INTERACTIONS DUE TO QUANTUM FLUCTUATIONS 



configuration n to be defined by (tib)x(t) = ti x (t) for any x x r G B, we introduce the useful 
notion of the restriction o>b of a quantum configuration u to B as to consist of ns and those 
quantum transitions from u that are contained in B, A x t C B (we suppose here that uj and 
B are such that no transition intersects both B and its complement; we do not define cjb in 
this case). 

Now the loops and and the quantum contours can be defined. First, the loops of a quantum 
configuration u are the triplets £ = (B,ub, g^)] B = Ax [ri,7^] G £?o(<*>) is the support of 
the loop £ and = n^Ti — 0) = tm.(t2 + 0), a restriction of a configuration g G G. (While 
the configuration g is not unique, its restriction to A is determined by the loop £ in a unique 
way.) We say that £ is immersed in g. Given a quantum configuration u>, we obtain a new 
configuration u by erasing all loops (B,UB,g^), i-e. for each £ we remove all the transitions 
in its support B and change the space-time configuration on B into g G G into which £ is 
immersed. Let us remark that B(6j) = Bi(u). Notice that, since we started our construction 
from (4.4), we have automatically diamA > 2Rq for a support ^4 x [ri,7^] of any loop £. 

Quantum contours of a configuration u will be constructed by extending pairs (B,u>b) 
with -B G B[(u) by including also the regions of nondominating states from G. Namely, 
summing over loops we will see that "loop free energy" favours the regions with dominating 
configurations from D C G. However, to recognize the influence of loops, we have to look on 
regions of size comparable to the size of loops. This motivates the following definitions with 
V(x) = {y G Z, u , \x — y\ < R} being an extension of original neighbourhood U(x). Thus, we 
enlarge the set E{n) of classically excited sites to E(n), with 

E(n) = {ieZ": n y(x .) / g v{x) for any g G G} 

and we introduce the set F(n) of softly excited sites by 

F(n) = {x G Z, u \ E(n) : nyi x \ ^ dyt x \ for any d G D}. 

Then, for a quantum configuration such that u = &, we define the new extended boundary 



B.(tf)= U ([E(n(T))UF(n(T))]xr)UU([U A V(x)] 

re [0,(3] per ' 1 = 1 



X T, 



and if u; ^ u, we set B c (u) = B c {u). Notice that B{Cj) C B e (u), since the first set 
is the union of classical excitations, quantum transitions and boxes; obviously the classical 
excitations and the quantum transitions also belong to B c (u), and the boxes being such that 
their diameter is smaller than R and they contain [/(x)-excited sites at each time, they are 
y(x)-excited. Decomposing B e (u) into connected components, we get our quantum contours, 
namely 7 = (B,ujb)- Notice that the configuration ujb contains actually also the information 
determining which dominant ground state lies outside B. We call the set B the support of 7, 
B = supp7, and introduce also its "truly excited part", the core, core 7 C supp7, by taking 

core 7 = supp7n(u re[0 , / 3] pcr (^(ri(-7")) x r) UU™^^, V(x)] X V Finally, notice that if 

the contour is not wrapped around the torus in its spatial direction, there exists a space-time 
configuration <jlP and we have B = B c (ui^). 

A set of quantum contours T = {71, . . . ,7*,} is called admissible if there exists a quantum 
configuration cj r G Q\° v which has T as set of quantum contours. Clearly, if it exists, it is 
unique under assumption that it contains no loops (u r = ui T ). We use P^ er to denote the set 
of all collections T of admissible quantum contours. 

Given T G V^ eT , a set of loops S = {£1, ■ ■ ■ is said admissible and compatible with T 
if there exists u; ru= which has H as set of loops and T as set of quantum contours (it is also 
unique whenever it exists). More explicitly, 
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• two loops £ = (B,UB,g\) and £' = (B',u B ,,g A ,) are compatible, £ ~ £', iff 5 U £>' is not 
connected; 

• using coreT = U 7g rcore7, a loop £ = (£?, w^, g^), with 5 = A x [n, T2], is compatible 
with r, £ ~ r, iff 

B U coreT is not connected, (4-5) 
gt A = n r A (T 1 -0)=n r A (T 2 + 0y, (4.6) 

• a collection of loops 3 = {£i, . . . , £^ } is admissible and compatible with T iff any two loops 
from H are compatible and each loop from H is compatible with T. 

We use 2?^ op (r) to denote the set of all admissible collections H that are compatible with T. 

The conditions of admissibility and compatibility above can be, for any given set of transi- 
tions {A\, . . . , A m }, formulated as a finite number of restrictions on corresponding transition 
times {ti, . . . , r m }. Given the restrictions on admissibility of T € T>^ cr , the restrictions on H to 
belong to V l ° op (V) fact orize. As a result, the partition function Z^ cr in ( f4.3|) can be rewritten 
in terms of integrations over T>^ er and J) l ° op (r) [the summation over T and S accompanied 
with the integration, a priori over the interval [0, /?], over times Tj of corresponding transitions, 



subjected to above formulated restrictions, c.f. (^2)]. Furthermore the contribution of V U 



factorizes as a contribution of T times a product of terms for £ € 5 [ BKU , DFF1]3 we get 



Z£ er = I dr / dEp^(u ru ^ 



pf 1 Jv l ° op (r) 



/ di>p-> r ) / d»n^)- 



'A 'CD £ £ S 

Here, using {(Aj, Ti), i = 1, . . . , m} to denote the quantum transitions of T U E, we put 



(4.7) 



p p Cr( ^ruH ) = j| (n ruH (ri _ 0)|FAjn ruH (ri + 0))exp |_ f a^^rus^^ (48) 

i=l ■'Ta 

where J B d(A, r) is the shorthand for J d r Z^A- J 4xrC-B ( use d here for 1? = Ta). Similarly for 
p pcr (uj r ). Further, the weight of a loop £ = (B^, u) B z, g A ) with the set of quantum transitions 
{{Ai, Ti),i = 1, . . . ,£} and the space-time configuration corresponding to Ugt, is 

z(0 = e*p{-J d(A,T)[^ A (n A (r)) - * A (g{)]} <4 1 V Al \n Ai ( n + 0)) x 

x (4 2 (r 2 - 0)| V A2 \n A2 (r 2 + 0)) . . . (n^fa - 0)| Va, |<4>. (4.9) 

Given T £ X^ er , ^ e secon d integral in (|4.7| ) is over the collections of the loops that interact 
only through a condition of non-intersection. This is the usual framework for applying the 
cluster expansion of polymers. The only technical difficulty is that the set of our loops is 
uncountable (the loops depend on continuous transition times), and thus we cannot simply 
quote the existing literature. Nevertheless, the needed extension is rather straightforward and 
often implicitly used. 



For spin or boson systems factorization is true simply because any two operators with disjoint supports 
commute. In the case of fermion systems there is an additional sign due to anticommutation relations between 
creation and annihilatio n opera tors, and factorization is no more obvious. That it indeed factorizes was nicely 
proved in Section 4.2 of jDFFlj]. 
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Given a collection C = (£1, . . . , £ n ) of loops, we define the truncated function 

$ T (C) = V(C) [J (4-10) 



n! 

«6C 



with 

V9 T (C) = ^ T (6,---,a 




if n = 1, 
- 1) if n > 2, 



where the sum is over all connected graphs Q of n vertices. Notice that <3? T (C) = whenever 
C is not a cluster, i.e. if the union of the supports of its loops is not connected. We use C\ and 
C\ to denote the set of all loops and clusters, respectively, and use J c dC as a shorthand for 
> i fctv ' ' Jc d£n, m ODV i° us meaning. Whenever T S V^ ei is fixed, we use C\(T) to 
denote the set of all loops compatible with T and write C S C\(T) whenever the cluster C con- 
tains only loops from C\(T). Again, Jc A (r) dC is a shorthand for ^ n > x Jc A tr) ' ' i*£ A (r) c ^ n ' 
Finally, we also need similar integrals conditioned by the time of the first transition encoun- 
tered in the loop £ or the cluster C. Namely, using C to denote the support of C, i.e. the 
union of the supports of the loops of C, and Iq = {ti(C), T2(C)} to denote its vertical 
projection^, Ic = {t G [0, /?] per ; Z" xrnC/f)}, we use Cjf' r) for the set of all clusters 
C £ Ca with the first transition time t\(C) = r, for which their first loop £i with support 
B\ = A\ x [ti(C),T2], contains the site x, Ai 9 x. Then f r (x,r) d£ and f r (x,r) dC are shorthands 

C A 

for the corresponding integrals with first transition time fixed — formally one replaces j d£i 
by J I \ A\ 3 x](5(ri(£i) — r)d£i. With this notation we can formulate the cluster expansion 
lemma. 

Lemma 4.1 (Cluster expansion). For any c G R, «i < (2Rq)^ i/ , ct2 < R' 2i/ Aq and 5 > 0, 
i/iere exists Eq > suc/t t/iai whenever \\V\\ ^ £o and T £ D P er ; we /iaue t/ie /oop cluster 
expansion, 



[ dHTTz(£) = exp{/ dC<£ T (C)l. (4.11) 

yj,ioop (r) Uc A (r) J 

Moreover, the weights of the clusters are exponentially decaying (uniformly in A and (5): 

f dCl[C3 (x,r)]|$ T (C)| 11 e (^Uog || v||)|A|+a 2 |B| ^ (4J2) 



and 



L 



dC|$ T (C)| TT e ( c - ailo sll v ll)l A l+ Q2 l B l < 5 (4.13) 

€6C 



for every (x, r) £T 



A- 



Proo f . One can follow an y st andard reference concerning cluster expansions for continuum systems, for ex ample 
[Bry . We are using here Pfi| whose formulation is closer to our purpose. Assuming that inequality (4.13) holds 
true, we have a finite bound 



d6...d^ T te,...,60ini z (£0K S0\A\. (4.14) 

£ A (r)" <=1 



Lemma 4.1 



then follows from Lemma 3.1 of [Pfij. Let us turn to the proof of the two inequalities. Let 

f/£\ _ U^| e C<:-"llosl|V||)|i4|+aa|B| 



3 Again, if ri > T2, the segment [n, T2] C [0, /3] pC r contains the point = f3. 
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Skipping the conditions ~ V, we define 



I n = n\j d£ 1 l[B 1 3(x,T)] + [ dfc] / d6...dC n | V T (ei,...,| n )|JI/(6) 



(4.15) 



(it does not depend on (x,r) £ Ta). The lemma will be completed once we shall have establis hed that 
In ^ n!(^5) n (assuming that 5 $J 1; otherwise, we show that I n ^ nl/2 n ). From Lemma 3.4 of Pfi|, we 
get 



W T (£,i, • • • , fn)| < ^ I! 1 t Bl U B J connecte< i] • 

T tree on n vertices e(i,j)€.T 



(4.16) 



Denoting di, . . . , d n the incidence numbers of vertices 1, . . . , n, we first proceed with the integration on the loops 
j 1 for which dj = 1; in the tree T, such j shares an edge only with one vertex i. The incompatibility between 
and £j, with £ = (Si, , <?^\ ), Bi — Ai x [r} , 7"2 ], and similarly for £ j, means that either Bj U [At x ] is 

connected, or [Aj x r^] U B t is connected. Hence, the bound for the integral over the £j that are incompatible 
with & is 



d£j I [Bj n Bi connected] /(£,-) < 2v\Ai\ / d^ I [Bj B (x, r)] /(£,) + 2v\B., 



s: 2i/(|Ai[ + a[B 



i|) C / dfcl [Bj 9 (s,t)]/&) + ~ [ u , d^m)). (4.17) 



(The constant a has been introduced in order to match with the conditions of the next lemma). Then 
/„ < n(2u) n ~ 1 



T tree of n vertices 



T \ d^l[B 1 B(x,r)}+ d^l/^OdArl+alBi 

L -'£ A JCY' T > 



n 

j=2 



dj-l 1 



d&I[Bj 3 (a;,T)]/(&)(|Aj| + a|Bj|) + - / ^ dfc/(&)( |Aj| + a|B, 



(4.18) 



Now summing over all trees, knowing that the number of trees with n vertices and incidence numbers di, ■ ■ ■ , d n 
is equal to 

(n-2)! 



(n-1)! 

(di - 1)! . . . (dn - 1)! " <*i!(da - 1)! • • . (dn - 1)! : 



we find a bound 

/„ < n\{2v) n -\l + a) 



dfl[B3(x,T)]/(0 



JA|+ Q |S| , j_ 



dC/(C)e 



A\+a\B\ 



(4.19) 



We conclude by using the following lemma which implies that the quantity between the brackets is small. 



□ 



Lemma 4.2. Let ct\ < (2Rq) v and 02 < R 2 ^Arj. For any c G R and 5 > 0, i/iere exists 
£o > suc/j t/iai whenever \\V\\ ^ £0 the following inequality holds true, 

I &tl\B3 (x,r)]|z(e)|e( c - ail °sll v ll)l A l +a2 l B l + / d£K£)| e^"" 1 l««IIV||)|A|+a 3 |B| < s 
Jc A Jct T) 

where (x,r) is any space-time site ofT\. 

Proof. Let us first consider the integral over £ such that its box contains a given space-time site. We denote by t\ 
the number of quantum transitions of £ at times bigger than r, and £2 the number of the other quantum transi- 
tions. The integral over £ can be done by summing over (£1 + £2) quantum transitions A{, . . . , A\ , A\ , . . . , Aj 2 , 

by summing over (£1 + £2) configurations , and by integrating over times T\ < • • ■ < t}^, T\ < ■ ■ ■ < t/ 2 . 

3 

Let us do the change of variables fx = t[ — t, t\ = t\ — r\ , ■ ■ ■ , t/ = t\ — t} 1 _x, and f 1 = r — rf, . . . , 
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ff = rf ! — Tg . Then we can write the following upper bound 



dp[_B3(x,r)]K6|e (c - ailog|ml)|A|+a2|B| < EE E / d n---cl<r! 



fa>l Ai,...Aj a „M,...,„W af£Gi! 

A connected 



4 



II II IK' 1^4- |n^ +1 )|e (c - ailog||v|l)l ^ 1 e -nZ A >cAl*A>(nJ)-* A ,( 9A ,)] e f)K"a 2 (4 2Q) 

8=1,2 j=l 

where (?a £ Ga is the configuration in which the loop £ is immersed (if the construction does not lead to a 
possible loop, we find a bound by picking any qa £ Ga)- Remark that we neglected a constraint on the sum 
over configurations, namely n^' 1 = n^' 1 . It is useful to note that the sums over £i,£2 and over the quantum 
transitions are finite, ot herw ise they cannot constitute a loop. 



Using the definition (2^5) of the norm of a quantum interaction, we have 

\{n A \V A \n A )\< W\\ W . 

Furthermore 



A'CA 

,\A\ 



as claimed in Property (2.3). Hence we have, since the number of configurations on A is bounded with S'~ 

, h r|| V -||l-a 1 (2flo)" n c(2R„)» 1 \A)\ 

/ d «i[fl 9 o B> r)]Koie^ ta -^'+-' fl '< e e n n [ " - R -, Ao R , a2 1 ■ 

Ca *1,*2>1 A},...A 2 e i=iaj=l 

A connected 

(4.21) 

This is a small quantity since the sums are finite, by taking ||V|| small enough. Now we turn to the second term, 
namely 



The proof is similar; we first sum over the number of transitions £, then over £ transitions A\,...Ai with 
A = VJiAi B x, A connected. Then we choose £ — 1 intermediate configurations. Finally, we integrate over £ — 1 



time intervals. The resulting equation looks very close to (4.20) and is small for the same reasons. 

□ 

Now, we single out the class of small clusters. Namely, a cluster is small if the sequence of its 
quantum transitions belongs to the list S. To be more precise, we have to specify the order of 
transitions: considering a cluster C = (£i, . . . , and using S(fW), £ = 1, . . . , Jfc, to denote the 

sequence of quantum transitions of the loop £^ = (5^,0)^,^ <S(£^) = S(B^\ w B ;<)), 
we take the sequence S(C) obtained by combining the sequences ),..., in this 

order. A cluster C is said to be small if S(C) € S, it is large otherwise. We use C^ ma11 to 
denote the set of all small clusters on the torus Ta- 

The local contribution to the energy at time r, when the system is in a state n^(r), is 
&A( n A( T ))- Similarly, we will introduce the local contribution of loops (and small clusters of 
loops) in the expansion of the partition function — the effective potential '^(n^r)). The 
latter is a local quantity in the sense that it depends on n only on the set A at time r. An 
explicit expression of ^ AgA) with g £ G is, in terms of small clusters, 

*5(<m) == - / dC, nTT^ 1 i c ~ 9a, a c =a,i c bo]. (4.22) 

7C| mal1 \*c\ 

Here, again, C is the support of C, Ac its horizontal projection onto Z", Ac = {x £ Tf\x x 
[0,/3] P er flC / 0}, and Ic its vertical projection, \Aq\ and \Ic\ their corresponding areas, and 
the condition C ~ gA means that each loop of C is immersed in the ground state g. Notice that 
"horizontal extension" of any small cluster is at most R: if C is a small cluster, diam (Ac) ^ R- 
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The definitions of Section |2.2| are now clear, once we identify the effective potential $ defined 
in (U) as the limit -> oo of ( p~22l) . Namely, 

^ = lim vL^. 

Our assumptions in Section 2.4 concern the limit (3 — > oo of the effective potential, but at 
non zero temperature we have to work with . To trace down the difference, we introduce 
= i$P _ Notice that ( [4.22] ) implies ^(n^) = whenever ^ GU or diamA < 2R . 

Recalling that if C C T\, C is the smallest box containing C, we introduce, for any cluster 
C G q ma11 , the function 

$T (c . r)= ^_iy drl [J C 9 r] (l [C ~ T] - I [n^(r) e G, C ,C ~ n^M]). (4.23) 

Here, the first indicator function in the parenthesis singles out the clusters whose each loop is 
compatible with T, while the second indicator concerns the clusters for which nF A (r) G Ga, 



if 



and each their loop is immersed in the configuration ti t a (t) (extended as a constant to all the 
time interval Ic)- Observing that <3? T (C;r) = whenever C ncoreT = 0, we split the integral 
over small clusters into its bulk part expressed in terms of the effective potential and boundary 
terms "decorating" the quantum contours from V. 

Lemma 4.3. For any fixed V G Da, one has 

[ dC<£ T (C) = -/ d(Ar)^K(r))-/ d(A, t)^J«(t))+ / dC<S> T (C;T). 

iq mall (r) Jj a Jt a JC| ma11 

The term 3> T (C;r) vanishes whenever C n coreT = 0. 
Proof. To get the equality of integrals, it is enough to rewrite 



L 



&C$ T {C)= dC$ T (C)I \C ~ rl (4.24) 

au (r) J cr M 



and 

- [ d(A,T)^ A (n r A (r))= [ dC$ T (C) rdrl[n r Ac (r)6G, c ,C~n^(r)]. (4.25) 
Jt a Jc^ ma11 Jo 

Moreover, whenever CflcoreF = 0, the configuration ti i Ac (t) belongs to Ga c , and it is constant, for all r £ Ic- 
Under these circumstances, the condition C ~ F is equivalent to C ~ ^a c (t") and the right hand side of (4.23) 
vanishes. 

□ 

Whenever V € T>\ is fixed, let W d (T) C Ta be the set of space-time sites in the state d, i.e. 

W d (T) = {(x,r) G T A : n T v{x) (r) = d v(x) }. 

Notice that 

T A = su PP ru u w d (ry, w d (r) n w d ,(r) = if d ^ d', 

and the set supp rnWd(r) is of measure zero (with respect to the measure d(x, r) on Ta). Let us 
recall that the equivalent potential T satisfies the equality ^AcA^( n ^) = SAcA(^( n ^0 + 
^ a(^a)) for any configuration n on the torus A and that we defined e(d) = YIabo Tj |a| 4 ' > f° r 
every d G D. 



Lemma 4.4. The partition function ( |4.7| ) can be rewritten as 

Z A per = f dr \{ e ^ r ^ n z( 7 ) e^ r ) 
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Here the weight 2(7) of a quantum contour 7 = (B,u>b) with the sequence of transitions 
(Ai, . . . , A m ) at times (n, ... ,r m ) is 

^(7)=nK 1 (^-0)l^ l |nX ! (r J + 0))exp{- / d(x, r)ej(n»)}. (4.26) 
i=i ^ 

T/te rest 7£(T) is giwen 6?/ 

K(F)= [ dC<J? T (C)-/ d(A,r)^(n^(r))+ / dC$ T (C ; r). (4.27) 

JC A (r)\c^ mall (r) JTa Jei mM 

Proof. Using the Lemmas and 4.3 to substitute in ( }4.7j ) the contribution of loops by the action of the effective 
potential, we get 

- m „ 

Zl CT = I per dr{n(nA 4 (r i -0)|Vk i |n^(r i + 0))}exp{- / d(A, T)(* A (r£(r)) + *a(»&(t)))} e TC(r) . 

(4.28) 

Replacing $ + ^ by the physically equivalent potential T, we get 

/ Tipcr dr{n(n r 4i (r l -0)|^ l |n^(r 1+ 0)>}e X p{- f ^ d(x, r)ej(n»)} ft *-'™ W '<™ e K < r > . 

(4.29) 



uppr 



We get our lemma by observing that the product over quantum transitions and the first exponential factorize 
with respect to the quantum contours, as it was the case for t he loop s (for fermions the sign is arising because 
of anticommutation relations also factorize; we again refer to [DFF1| for the proof). 

□ 

Our goal is to obtain a classical lattice system in v + 1 dimensions. Thus we introduce a 
discretization of the continuous time direction, by choosing suitable parameters (3 > and 
N E N with = Setting La to be the (y + l)-dimensional discrete torus La = 

A x {0, 1, . . . , N — l} per — let us recall that A has periodic boundary conditions in all spatial 
directions — and using C(x,t) C R^ +1 to denote, for any (x,t) £ La, the cell centered in 

(x, ^t) with vertical length /3/A, we have Ta = ^( x ,t)£h A C(x,t). 

For any M C La, we set C(M) to be the union of all cells centered at sites of M, C{M) = 
U( x ,t)eMC(x,t) C Ta- Conversely, if B C Ta, we take M(B) C La to be the smallest set such 
that C(M(B)) D B. Given a connected! 12 ] set M C La and a collection of quantum contours 
T G V^ er , we define 



ip(M; r) = / dCI [M(C) = M] $ T (C) + 

•/c A (r)\c? mall (r) 

f/ dCl[M(C) = M,C , £C7(sup P r)]$ T (C;r)- f d(A,r)^5(n£(r)) (4.30) 

JC? ma11 JM(Axt)=M 



and 



K(T)= [ dCI [C C C7(su PP r)]$ T (C;r). (4.31) 

J £ small 

We have separated the contributions of the small clusters inside C(suppr), because they are 
not necessarily a small quantity, and it is impossible to expand them. On the contrary, <p(M; T) 



11 Remark the difference from [BKU]; here the vertical length of a unit cell /3/A depends on ||V"||, since so 
does the quantum Peierls constant A. 

12 Connectedness in La is meant in standard way via nearest neighbours. 
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is small, and hence it is natural to write 



£ II (e^ (M;r) -l), (4.32) 
M MeM 

with the sum running over all collections A4 of connected subsets of La- 

Let suppA4 = UmemM. Given a set of quantum contours T 6 T>^ T and a collection Ai, 
we introduce contours on La by decomposing the set M(suppT) U suppAl into connected 
components [notice that if (x,t) ^ M(suppT) UsuppA4, then C(x,t) C ^deDWd(T)]. Namely, 
a contour Y is a pair (suppY, ay) where suppY C La is a (non-empty) connected subset of 
La, and ay is a labeling of connected components F of dC(suppY), ay(F) = 1, . . . ,r. We 
write |Y| for the length (area) of the contour Y, i.e. the number of sites in suppY. A set 
of contours 3^ = {Yi, . . . ,Yk} is admissible if the contours are mutually disjoint and if the 
labeling is constant on the boundary of each connected component of Ta \ Uy e yC(suppY). 
Finally, given an admissible set of contours y, we define Wd(y) to be the union of all connected 
components M of La \ Uygysupp Y such that C(M) has label d on its boundary. 

Consider now any configuration u) € QF r yielding, together with a collection Ai, a fixed set 
of contours y. Summing over all such configurations u and collections Ai, we get the weight 
to be attributed to the set y. Let r w be the collection of quantum contours corresponding 
to u>, UygysuppY = M(suppT a ') UsuppA4. Given that the configurations u are necessarily 
constant with no transition on Ta \ C(M(suppT w ) U suppTW), we easily see that the weight 
factor splits into product of weight factors of single contours Y € y. Namely, for the weight 3 
of a contour Y we get the expression 



3(F) = / dr FT z(i) TT e -<d)\w d (r)nc(su PP Y)\ e 7e(r) 
JvriY) 7 er den 

1 [M(supp T) U supp M = supp Y) Yl ( e v(M;r) - l) , (4.33) 
M MeM 

where T>? eT (Y) is the set of collections V of quantum contours compatible with Y, T E T>^ r {Y) 
if suppT C supp Y and the labels on the boundary of suppT match with labels of Y . Thus, we 
can finally rewrite the partition function in a form that agrees with the standard Pirogov-Sinai 
setting, namely 

Z F = E II e-l^lWI I] 3(n (4-34) 
y deD Yey 

with the sum being over all admissible sets of contours on La- 

In the next section we will evaluate the decay rate of contours weights in a preparation to 



apply, in Section y, the Pirogov-Sinai theory to prove Theorems BTL 2.2, and 2.3. 



5. Exponential decay of the weight of the contours 

In this section we show that the weight 3 has exponential decay with respect to the length 
of the contours. We begin by a lemma proving that the contribution of A4 is small, that we 



shall use in Lemma 5.2 below for the bound of 3. 

Lemma 5.1. Under the Assumptions for any c < 00 there exist constants (3q, (3q < 00, 
and £q > such that for any (3 ^ (3q, (3q ^ f3 < 2(3q, and ||V|| ^ Eq, one has 



E 



M3[x,t) 



for any contour Y and any set of quantum contours T E P^ er (Y). 
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Proof. We show that 



This implies that |<^(Af;r)| ^ 1 and c onsequently Lemma 5.1 



Y, \f(M;F)\e c]Ml < 1. 

M3(i,t) 

holds - 



with a slightly smaller constant c. 



Let us consider separately, in ( 4.3C ), the three terms on the right hand side: (a) the integral over big clusters, 
(b) the integral over small clusters, and (c) the expression involving t/j 13 . 
(a) Big clusters. Our aim is to estimate 



J= £ « 

M3(i,() 



c\M\ 



c A (r)\c?»=-»(r) 



AC I [M(C) = M] |$ T (C)|. 



Since M(C) — M and M 9 (x, t), the cell C(x, t) either intersects a quantum transition of C, or it is contained 
in a box B belonging to a loop of C (both possibilities may occur at the same time). In the first case we start 
the integral over clusters by choosing the time for the first quantum transition, which yields a factor /3/A. In 
the second case we simply integrate over all loops containing the given site. In the same time, given a cluster 
C = (£i, . . .,£„), & = (Bi.wg^p^.) and B t = Ai x [t^\ t 2 W ], the condition M(C) = M implies that 



A 

^jl^l + ll^l] > |A/|. 



Using it to bound \M\, we get the estimate 

c|A|+c4|B| 



dC|$ T (C)| J] e 

Sec 



+ 



dCl[C3(x,r)]\<£ T (C)\ fj c 



c|A|+e# |fl| 



(5.1) 



(5.2) 



Taking, in Lemma 4.1, the constant c as above as well as a\ = j(2i?o) Q2 = cA//3, 5 = 1, and choosing 
the corresponding £o(c, ai, aa, 5), we can bound the second term of (|5.2|), for any ||V|| ^ £o, with the help of 
(1.12) once (3 is chosen large enough to satisfy 



A A 



(5.3) 



To estimate the first term of (5.2), we first consider the contribution of those clusters for which 



£ 

A 



{2R )-"\A\ 



Applying it together with (5.3) we can directly use the bound (4.13). 
Thus it remains to estimate the contribution of those terms for which 

2 - E w < los(A//5) 



(2i?o) 



fee 



log||V|| 



Let us first fix (3 and eo ^ eo(c, ai, «2, 5) with the constants c, ai, 02, and 5 as above, so that 



and, in the same time, 



i- > ^-B?» 
so A 



a , fc-2fc'(2fl )- 



(5.4) 



(5.5) 



(5.6) 



for a suitable large k' (we also assume that £0 ^ 1). Here k is the constant that appears in Assumption ^, 
A(||V||) > \\V\\ ■ Observing further that A(||V||) can be taken to increase with ||V|| (one can always consider 
a weaker lower bound A when taking smaller ||V||), we conclude that (5.S), as well as the condition 



,log(A//3) 



< fc'. 



log || V|| 

are satisfied for every ||V|| ^ eo- Thus, it suffices to find an upper bound to 



J'-I 
A 



(5.7) 



dC|* T (C)|lE \A\<k'}. 
eec 

The main problem in estimating this term stems from the factor 1/A that may be large if ||V|| is small. Thus, 
to have a bound valid for all small ||V||, some terms, coming from the integral, that would suppress this factor 
must be displayed. 
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The condition X^ec \A\ < k' wn l be used several times by applying its obvious consequences: (i) the number 
of loops in C is smaller than k' , (ii) the number of transitions for each loop is smaller than k' , (iii) each transition 
A is such that \A\ < k' , and (iv) the distance between each transition and x is smaller than k' . 

Furthermore, we use Assumption ^ to bound the contribution of the transitions of C ; recalling the definition 
(fLSl) of the weight of £, we have, for any large C, 



II < 6 i(l|V||)A II ex P{" / d(A,T)[* A (nl(T)) - <Mffi(r))]} 



fee 



(5.8) 



5ec 



In the last inequality we used Assumption |^ in the form of the bound ( 2^3 ) as well as the lower bound |ts — ri | = 
jjj ^ for the support B — Ax [n.,75] of the loop £. 

For any £ 6 C = (£1, . . . £n), let r be the time at which the first transition in C occurs (we assume that it 
happens for the "first" loop £1) and be such that r + is the time at which the first transition in £ occurs 
(t ?1 = 0). Referring to the condition (i) on the number of loops in C, we get the inequality 

£ |r«Kfc'£|B|, 



and thus also 



1 TT e iy r' 1 

e, 



r|T«| 



Integrating now over the time of the first transition for each £ 6 C, £ ^ fi, and taking into account that 



|v A (ft,---.WI<n" 



we get 



?2k'R 2 



d£e~ 



'A |S| 



■ 



(5.9) 



Here the constraint I [£j : fc'l means that the loop £j satisfies the conditions (ii)-(iv) above. We have then a 
finite number of finite terms, the contribution of which is bounded by a fixed number K < 00 (depending on eo, 
$, and k'). Thus J' ^ /9&i(||V||)if which can be made small by taking ||V|| sufficiently small. 

(b) Small clusters. Let us first notice that $ T (C;r) < $ T (C)j, and since M{C) = M, inequality ( |5.l|) is 
valid. Moreover C must contain at least one of the two boundary points (y,t-^ ± ^v) of some cell C(y,t) for 
which dist (x, y) ^ R. Indeed, given that C is s mall a nd in the same time CTlcore r j= (c.f. Lemma (i|), this is 
the only way to satisfy also C C(suppF) [c.f. (4.30)]. Thus it suffices to use again (]4. 12|) and (5.3) to estimate 



{2ny 



ACl\C3{x,r)]\^{C)\\\ e cW+c 7 W . 

In small 

' c a tec 
(c) Bound for tp 13 . Finally, we estimate the expression involving i/A We first observe that 

e° ? l<(9A)Kl 
for any AcZ" and with a = |f? _2l/ Ao, Indeed, 



(5.10) 



Wa{9a)\ 



all 



¥ 3 A {g A )-H! A {g A )\ 



a/3 



/ ACl [C ~ g A , A c = A, I c 9 0, C C A x [0, p] pcr , |/ c | = P] ^rr=r + 

J c small \Io\ 

+ [ ACI\C ~ g A ,A c = A, la 3 0,C C A x [-00, 00], llcl > /3l 

Jcjmall 



$ T (C) 



(5.11) 

The first integral above corresponds to clusters wrapped around the torus in vertical direction, while the second 
one assumes integration over all clusters in A x [—00, 00]. For any C above it is \Ic\ zZ P and thus 

B af> ^ -Q e a|fl| 

tec 

Observing now that every cluster in both integrals necessarily contains in its support at least one of the points 
(x,0), x £ A, and using the fact that diamA ^ R, we can bound the first integral by 

dci[c 3 (x,o)]i$ T (c)in C ' 



,\B\ 
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which can be directly evaluated b y (|4.12j ). The same bound can be actually used also for the second integral, 
once we realize that the estimate ( |4.12[ ) is uniform in [3. 

Using now the fact that tf>^ = if diamyl ^ R, the condition M(A x {r}) = M implies that M has less than 
R" sites, hence e c ' A/ ^ e cR . Furthermore, referring to ( [5.10[ ), we have 

d(,4,r)j^(-)|l[M(Ax{r}) = M]e c|M| <: !L e -h R ' 2 " ^+<= R \ (5.12) 

which can be made small for j3 sufficiently large and concludes thus the proof of the lemma. 

□ 



Using Lemma 5A and introducing eo = min^eD e(d), we can estimate the weight 3 of the 
contours in the discrete space of cells. 

Lemma 5.2. Under the Assumptions for any c < 00, there exist (3q,$o < 00 and eo > 
such that for any (3 ^ (3q, /3q ^ f3 < 2(3q, and ||V|| ^ eo, one has 

| 3 (Y)| <; e ~l e °l y l e" c l y l 

for any contour Y . 

Proof. For a given V (such that M(supp F) C supp Y) with transitions {A±, . . . , A m } at times {n , . . . , r m }, we 
define A(T) = U™ 1 L> xeA , [V(x) x n], A = M(A(T)), and £ C suppF \ A to be the set of sites (x,t) such that 
n v(x)( T ) 4- Dy(x) for some (x,r) £ C(x,t). The latter can be split into two disjoint subsets, £ — £ COTC \j £ Bott ; 
with (x,t) £ ^c 010 whenever (1^1(7) ^ f° r some (a;,r) £ C(x,t). The condition M(suppF) UsuppA'f = 

suppY" in (4.33) implies the inequality 

e c\Y\ ^ e c(2Rr\A(r)\ e c\S\ -Q e c\M\ _ 

From definitions ( |4.33[ ) of a(Y) and ( |4.26| ) of 2(7), and using Assumption ^| we have 

e C ' y| 13(^)1 < E e"^ eo|SUPPrX " 41 V V e -(/3-c)|£\£-™| e _(|*p.(2K)-" - c )|£°°™ I x 

^icsuppy £csu PP y\^i£ core c£ 



|e c(2H)-|^l x 



r dri[M(A(r)) = AM(corer) =f coro ] ni«(^-0)|V Al |n^(ri + 0))|. 

i=l 

xexp{- f d(a;,r) e J(n r (r))}el^ r »l E [] | e^ M « - l| . (5.13) 

•'C(A') ' Al.supp MCsuppY MeM 

All elements in M are different, because it is so in the expansion ( |4.32| ). Therefore we have 

n l 6 ^"^ "l|e c|M| < E M E l e ^ M:r) -l|e c|M| ]" 

A-l,supp AlCsupp Y M<=M. n ^ ' MCsuppF 

< E ^[l Y l E |e^ M;r) -l|e c|M| l n (5.14) 

and us ing L emma |5.l| this may be bounded by e' y . 

In (4.31) clusters are small, and they must contain a space-time site (x,t) such that there exists x' with 
(x 1 , t) £ core Y and dist (x, x') < R. So we have the bound 



\H{F)\ < (27?)" I core T I / dC I [C 9 (x, r)] |<E> T (C)| , 

JC small 



since |$ T (C;r)| ^ |<E> T (C)|. Taking now, in Lemma 4.1, the constants c = a± — 02 = and 5 = 47^^" > an< ^ 
choosing the corresponding £0, we apply (4.12) to get, for any ||V|| ^ eo, the bound 

\TZ(r)\ < ^£(27?)-"|corer| < (2_R)-"|£ core |. 



Assuming f3 > c and f 4p ^ (2R) V c [c.f. Q], we bound e-^- c >l £ \ £0ar ° I e"^^ 2 ^ "- c )l£ coro I <; 1 
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Inserting these estimates into ( |5,13| ) , we get 

e e|V| | 3 (Y)|< e-^ eo|y[ e |y| ]T / drl [ M (A(r)) = A] ]J\ {n T A .(n - 0)| V Ai \n%{n +0)> | x 

-4csu PP y • /t 'a ct i=i 

x e <WI^I e r /" d (x,r)[ e J(n r (r))- eo ]). (5.15) 

1 JC(A) ' 

To estimate the above expression, we will split the "transition part" of the considered quantum contours into 
connected components, to be called fragments, and deal with them separately. Even though the weight of a 
quantum contour cannot be partitioned into the corresponding fragments, we will get an upper bound combined 
from fragment bounds. Consider thus the set 

A{T) = corernC(A(r)) 

and the fragments = (Bj,ws ( ) on the components Bt of A(F), ACT) — Uf =1 Bi, u>Bi is the restriction of o; r 
onto Bi. 

From Assumption W, we have 



d(x,r)[eJ(n r (r))- eo ] > \{2R)-» A 1*1 



Let us introduce a bound for the contribution of a fragment £ with transitions Aj, j — 1, . . . ,k, 

/. 

m = e -4(«D-"A |B| -q |( n C.( Tl _o)|^. |„C .( Tl + )>|eW^I . 
3=1 

Then, integrating over the set J-c(A) of all fragments in C(A), we get 

e c l y l \i(Y)\ < e-^ e °l y l e ly| £ 3^"^ £ A ( / W))* (5-16) 



.ACsuppy n > • / - A C(A) 



Anticipating the bound \ T d£f(£) ^ |.A|, we immediately get the claim, 

e c|y| | 3 (Y)| < e-^ eo|yl e 3|y| , 

with a slight change of constant c — > c — 3. 

^4 bound on the integral of fragments. Let us first consider short fragments £ = (B,u:b) satisfying the 
condition 

3=1 

The integral over the time of occurrence of the first transition yields the factor $/A. Notice that £ is not a loop. 
This follows from the construction of quantum contours and the fact that B is a connected component of A(T), 
where every transition is taken together with its /^-neighbourhood. Thus, either its sequence of transitions does 
not belong to <S, or the starting configuration does not coincide wi th th e ending configuration. In the first case 
we use Assumption in the second case Assumption ^, and since (5.17) means that the sum over transitions is 
bounded, we can write 

f dC8(C)<§l-4- (5-18) 
J:F c\% 

Finally, we estimate the integral over £'s that are not short. We have 

/ dC«(C)<Hl4 /, , dCi(C). (5-19) 

Here is the set of all fragments £ whose first quantum transition (Ai,n) is such that x £ A\ and r = Ti. 

Whenever C, is not short, we have 

i^niivir^' 1 . 

p j=l 

Thus, defining 



z'(0 = e-iwr^ow -Qj|| y ||i e c ( 2Kr+ipi 5 (5 20) 
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we find the bound 



\A\ / dCi'(C)- 



Here, slightly overestimating, we take for T(x, r) the set of all fragments containing a quantum transition (A, r) 
with x £ A. 

The support B of a fragment £ = (B,ujb) £ J-(x,t), is a finite union of vertical segments (i.e. sets of the 
form {y} x [n, T2] C Ta) and k horizontal quantum transitions A\, . . . , A^. 
We will finish the proof by proving by induction the bound 

/ dCi'(C) < 1 (5.21) 

with J-(x, t; k) denoting the set of fragments from J-(x, r) with at most k quantum transitions. 

Consider thus a fragment ( with k horizontal quantum transitions connected by vertical segments. Let (A, r) 
be the transition containing the point (x,t) and let {A\,t + n), . . . , (Ae,r + re) be the transitions that are 
connected by (one or several) vertical segments of the respective lengths |ti|, . . . , \re\ with the transition (A, r). 
If we remove all those segments, the fragment £ will split into the "naked" transition (A, r) and additional £ < I 
fragments Cl, . . . , Cl> such that each fragment Q, j — I, ...,£, belongs to J-(yj , r + Tj ; k — 1) with yj 6 A. Taking 
into account that the number of configurations (determining the possible vertical segments attached to A) above 
and below A is bounded by S 2 ^ A ^ and that the number of possibilities to choose the points yj is bounded by 
\A\ e , we get 



/ dCz'(C) < 

A.dist (A,x)<R 1=1 J J 3=1 J ^(Vj,r + Tj;k-l) 

< E [|[V||*^e c(2 ^ +2 ] |A| e a < a *>" /Ao <1 ( 5 - 22 ) 

A.dist (A,x)<R 

once ||V|| is sufficiently small. □ 

In the application of Pirogov-Sinai theory we shall also need a bound on derivatives of the 
weight of contours. 



Lemma 5.3. Under the Assumptions for any c < 00, there exist constants a,/3o,f3o < 00 

I - 

dm 



and £q > such that if (3 ^ [3q, Po ^ (3 < 2/3q, and [|V|| + X)i=i ll^^ll ^ £ 0; one /10s 



for any contour Y . 



Proof. From the definition (4.33) of 3, one has 

l£Kni<wni{EI^ 

dr [*( 7 )| n e -e"Wl^nc(suppr)| e !* ( r)| 
^ I [M(sup P r) UsuppM = suppF] £ J e v(M;r) A^( M; r)| JJ | e v(M ' ;r) - l|. (5.23) 

The bound for |g^_z(r)| is standard, see [BKU], and \^-e ti (d) \ is assumed to be bounded in Assumption |^. For 
the other terms we have to control clusters of loops. Since we have exponential decay for z(£) with any strength 
(by taking f3 large and ||V|| small), we have the same for g^rz(£) (by taking (3 larger and ||V|| smaller). The 
integrals over C can be estimated as before, the only effect of the derivative being an extra factor n (when the 
clusters have n loops). □ 
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6. Expectation values of local observables and construction of pure states 

So far we have obtained an expression (|4.34j) for the partition function Z? ei of the quantum 
model on torus A in terms of that of a classical lattice contour model with the weights of the 
contours showing an exponential decay with respect to their length. Using the same weights 
3(Y), we can also introduce the partition functions ^a(L) W1 ^h the torus A replaced by a 
hypercube A(L) and with fixed boundary conditions d. Namely, we take simply the sum only 
over those collections y of contours whose external contours are labeled by d and are not 
close to the boundary^ . Notice, however, that here we are defining ^(L) directly in terms 
of the classical contour model, without ensuring existence of corresponding partition function 
directly for the original model. We will use these partition functions only as a tool for proving 
our Theorems that are stated directly in terms of quantum models. 

To be more precise, we can extend the definition even more and consider, instead of the torus 
A, any finite set V C L = 7L V x {0, 1, . . . ,N — l} per . There is a class of contours that can be 
viewed as having their support contained in V C L. For any such contour Y we introduce its 
interior Int Y as the union of all finite components of L \ supp Y and Int jY as the union of all 
components of Int Y whose boundary is labeled by d. Recalling that we assumed v ^ 2, we note 
that the set L \ (suppY U Int Y) is a connected set, implying that the label ay(-) is constant 
on the boundary of the set V(Y) = supp Y U Int Y. We say that Y is a d-contour, if ay = d 
on this boundary. Two contours Y and Y' are called mutually external if V(Y) n V(Y') = 0. 
Given an admissible set y of contours, we say that Y E y is an external contour in y, if 
supp Y n V(Y') = for all Y' <=y,Y' + Y. The sets y contributing to Z^ are such that all 
their external contours are d-contours and dist (Y, dV) > 1 for every Y € y. 

In this way we found ourselves exactly in the setting of standard Pirogov-Sinai theory, or 
rather, the reformulation for "thin slab" (cylinder L of fixed temporal size N) as presented 



in Sections 5-7 and Appendix of [BKU]. In particular, for sufficiently large (3 and sufficiently 



small ||^||+^i=i llzfr^ll' tli ere exist functions /^(d), metastable free energies, such that the 
condition Kef /3 ' fl (d) = /o, with f = defined by /o = mhv e _D Re /^'**(c£')j characterizes 
the existence of pure stable phase d. Namely, as will be shown next, a pure stable phase {•)% 
exists and is close to the pure ground state \d). 

There is one subtlety in the definition of f^ ,fl (d). Namely, after choosing a suitable (3q, given 
(3, there exist several pairs (j3,N) such that € (/?o,2/3o) and Nj3 = (3. To be specific, we 
may agree to choose among them that one with maximal N . The function fP'^{d) is then 
uniquely defined for each (3 > (5q. Notice, however, that while increasing /?, we pass, at the 
particular value (3m = N(3q, from discretization of temporal size N to N + 1. As a result, the 
function f^ ,fl (d) might be discontinuous at (3n with [3 = co being an accumulation point of 
such discontinuities. Nevertheless, these discontinuities are harmless. They can appear only 
when Re f^'^(d) > fo and do not change anything in the following argument. 

Before we come to the construction of pure stable phases, notice that the first claim of 
Theorem |2.2| (equality of /o with the limiting free energy) is now a direct consequence of the 
bound 



(c.f. |BKUj , (7.14)). Here Q = {cZ; Re /^(d) = / }. 



The expectation value of a local observable T is defined as 



/ - "'j^L, «->) 



13 In the terminology of Pirogov-Sinai theory we rather mean diluted partition functions — see the more 
precise definition below. 
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In Section |] we have obtained a contour expression for Z? ei = Tr e~P H l . We retrace here 
the same steps for Z^ er (T) := TrTe - ^ . Duhamel expansion (4.1) for ZF eT (T) leads to 
an equation analogous to ( [L2]) , 



Zr(T)=Yl E E / d7i...dr ra (ni|T|ni> 

n ™ A A J0<TT<...<T m <3 



C 



m > n r .. n m Ai,...,A m •^0<T 1 <...<r ra < / 3 
AiCA 

--< )PCr «) (ni| ^ |ni) e -(---)< )PC >i) . . . <„£| V Am |n A > e^-K^K) . (6.3) 



Configurations n\ and match on A \ supp T (supp T C A is a finite set due to the locality 
of T), but may differ on suppT if T is an operator with non zero off-diagonal terms. Let 
QF> cr (T) k e ^ e ge ^. q £ q Uan t um configurations with n\(r) that is constant except possibly at 
U™i(A x Tj) U (suppT x 0). Then 

Zr r (T)= / du TA (n° A \T\ni}p^(u TA ). (6.4) 

We identify loops with the same iteration scheme as in Section ||, starting with the set 
B(°\uj) U(suppT x 0) instead of B(°\u>) only. This leads to the set B T {u). Removing the 
loops, we define B^(uj), whose connected components form quantum contours. There is one 
special quantum contour, namely that which contains suppT x 0. Let us denote it by 7 T and 
define its weight [see ( |4.26| )1 

m 

- T (7 T ) = KT ppT (-0)| T |< T ppT (+0)) n<nj(7, - 0)| V Ai \ntM + 0)> 



i=i 



exp{-y d(x,r)ej(n^(r))}. (6.5) 



Let T T = {7 T , 71, • • • , 7fc} be an admissible set of quantum contours, defining a quantum 
configuration u r £ Q^ T (T). Then we have an expression similar to that of Lemma [4.4j , 



Zr(T)= [ dr T TT e-l^( rT )l £ M/( 7 T ) TT ^(7)e^ (rT) , (6.6) 

T>P CT (T\ 

A l J J dGD 7 grT\{ 7 T} 

with ^(r T ) as in (^2^ ) with T replaced by T T . 

Next step is to discretize the lattice, to expand e^( r ' , and if Y T is the contour that contains 
suppT x C L A , to define i T {Y T ) [see (gf|)]: 

3 T (o = / pr dr^ r ( 7 T ) n -(7) n ° 

a ( yT ) 7er T \{ 7 T } der» 

^ I [M(supp T T U supp M = supp Y T ] ]"[ ( e ^ (M;rT) - l) . (6.7) 
M MeM 

We also need a bound for 3 T (y T ). It is clear that the situation is the same as for Lemmas 



-e(d)\w d (r T )nc( S up P Y T )\ e ii(r T ) 



T T 

5.1 and |5.2| , except for a factor (n2 uppT (~ 0)| T l n J U p P r(+0)) that is bounded by ||T||. We can 
thus summarize: 

Lemma 6.1. Under the Assumptions for any c < oo, there exist (3o>0o < °°> an d £ o > 
such that if ^ /?o, /?o ^ /3 < 2/?o and ||V"|| ^ £o 3 we have 

zrm= e n e " |e(d)|w<i(yT)l 3 T (^ T ) n *(n ( 6 - 8 ) 
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for every local observable T , with 
for any contour Y T . 



\l T {Y T )\ < ||T||e c|suppT| e -^ eo|yT| e" c|yT| 



In a similar manner as at the beginning of this section, we can introduce Zy(T) for any VcL 
by restricting ourselves in the sum ( |6,8| ) to the collections y T whose all external contours are 
(i-contours and dist (Y, dV) > 1 for every Y G y T . Thus we can define the expectation value 

(T) d y = 4P (6-9) 

for any V C L and, in particular, the expectation (T)^^ for a hypercube A(L). 

Again, this is exactly the setting discussed in detail in [BKU]. We can use directly the 
corresponding results (c.f. |BKU ], Lemma 6.1) to prove first that the limiting state (•)« exists. 
Further, retracing the proof of Theorem 2.2 in JBKU ] we prove that the limit 

TrTe -/3f/P cr 

K/i" Tr e~@ H T 
exists for every local T (proving thus Theorem ^]l|). Moreover, 



mr = . h k.. ~ ( 6 - 10 ) 



where, again, Q denotes the set of stable phases, Q = {d; Re /^(d) = /o}- Thus we proved 
the claim d) of Theorem |2.2j . 

Also the assertion c) follows in standard manner from contour representation employing di- 
rectly the exponential decay of contour activities and corresponding cluster expansion (c.f. | [BKU |, 
(2.27)). 

Before passing to the proof of b), we shall verify that (-)'g is actually a pure stable state 
according to our definition, i.e. a limit of unperturbable states^. To this end, let us first 
discuss how metastable free energies f" ,fl (d) change with [i. The standard construction yields 
fp>v(d) in the form of a sum e^(d) + s^ ,fM (d), where s l3 ' fJ- (d) is the free energy of "truncated" 
contour model K' d (Y) (see [ BKUfl , (5.13) and (5.6)) constructed from labeled contour model 
(4.34), which is under control by cluster expansions. As a result, we have bounds of the form 

0(e~ ,9 + || V|| + YH=i ) on |s^' M (cf)| as well as on the derivatives with respect to fi. Hence, 
in view of Assumption 7, the leading behaviour is yielded by e^(d). 

Starting thus from a given potential with = {d € D; Re f^'^{d) = ftf}, one can easily 
add to 3?^ a suitable "external field" that favours a chosen d E Q^. For example, one can take 

^ X a (n) = ^ A (n)+a6 d A (n) 

with 5 A defined by taking 5 A (n) = for n A = d A and 5 A (n) = 1 otherwisd^- Now, since 
^da ^ * s bounded from below by a positive constant (while ^Q^f - = for d' ^ d), for any 
a > the only stable phase is d, Re f^'^' a (d) = fQ ,fl ' a = min^gz) Re/^ ,At,a ((i'), and, in the 
same time, Ref^' a (d') > ff^ a for d' + d. Thus, Q^ a = {d} and (-)%^ a = OJ^V This 
state is unperturbable — when adding any small perturbation, metastable free energies will 
change only a little and that one corresponding to the state d will still be the only one attaining 



14 Recall that, up to now, the state (•)« is defined only in terms of the contour representation (see (6.9), (6.< 
and (4.34)), and the only proven connection with a state of original quantum model is the equality (6.11). 
15 Actually, we can restrict 5% only to a particular type of sets A — for example all hypercubes of side R. 
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the minimum. The fact that in the limit of vanishing perturbation we recover as well 

as the fact that 

lim {-)% cv = lim (•)!,„= (■)£„, 

follows by inspecting the contour representations of the corresponding expectations and ob- 
serving that it can be expressed in terms of converging cluster expansions whose terms depend 
smoothly on a as well as on the additional perturbation. 

To prove, finally, the claim b) of Theorem ^L^, it suffices to show that it is valid for ( - )^^ Q = 
i^p^a for every a > 0. Abbreviating {-)^ a = (-} per and H^ a,per = H^ ei , we first notice 
that the expectation value of the projector onto the configuration d on suppT, Pg UppT ■= 

|4up P t)(4u PP t| , is close to 1, since its complement ((1 - Ps UppT )) per = (Q 1 ~ P su PP T)) d is 
related to the presence of a contour intersecting or surrounding suppT (loops intersecting 
suppT x {0} are considered here as part of quantum contours), whose weight is small. More 
precisely, for any 5 > we have 

((]l-^uppT)) PCr ^^UppT|, 



whenever ||V[| is small enough and (3 large enough. Furthermore, 
^pcr _ er T; r (-P S u ppr TP s ^ ppT e fSH ^ ) + 

+ Tr ((1 - Pi ppT )TPl ppT e" W ) + Tr (T(l - P s d uppT ) e" W )] (6.12) 

and 

Tr(P s t ppT TP s rf uppT e-^A Cr ) = (d A \ T |d A )Tr (P s d nppT e"^ ) 

= (d A \ T \d A ) [Tr ( e" W ) - Tr ((1 - Pi ppT ) e~ W )] , 

(6.13) 

so that we have 

\{T)Z ei -(d A \T\d A }\< 

< Kd A | t \d A )\((i - pi PpT ))r + |p- ^u PP T)^uppT)r r i + \<m - p^^n 

(6.14) 

The mapping (T, T") i— > (T^T') per , with any two local operators T,T', is a scalar product; 
therefore the Schwarz inequality yields 



| (T)r - (d A \ T \d A ) | < | (d A \ T \d A ) | ((1 - P s d uppT ))r 

+ (((i - ^u PP T)) A pcr ) 1 ( i(pi PP T^TPi PPT )r] h + [(rtr)n *) 



sS IITI 



<(1- ^uppt)^ 1 ' + 2(((1- P s t PP T))D 1/2 < ||r|||suppT|(5 + 2^). (6.15) 



The proof of the remaining Theorem 2.3 is a standard application of the implicit function 
theorem. Thus, for example, the point /xo of maximal coexistence, Re (<i) = Re/^'^ (d') 
for every pair d, d! G D, can be viewed as the solution of the vector equation f(flo) = 0, 
with /(//) = (Ref^(di) - Re/^(4))[=i- Now, / = e + s, e(ji) = (e"(di) - e^))^ 1 , 
s(/x) = (Re s^ ,fl (di) — Re s^'^{d r )) r i z\, with ||s|| as well as ||J^|| bounded by a small constant 

once ||V|| + Ya=i |zr i s sufficiently small (3 is sufficiently large. The existence of a unique 
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solution hq G U then follows once we notice the existence of the solution /j,q <G U of the equation 
e(/xo) = (equivalent with e tM> (d) = e tJ, °(d'), d,d' € D) and the fact that the mapping 

with A~ l the matrix inverse to (^), is a contraction. To this end it is enough just to recall 
Assumption 7 and the bounds on s l3 ^{d), d € D, and its derivatives. 
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